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A  Systems  Programmer's  Guide  For  Implementing  OMNI TAB  II 
Sally  T.  Peavy,  Ruth  N.  Vainer,  and  Shirley  G.  Bremer 


OMNITAB  II  is  a  general-purpose  program  which  permits  direct  use  of 
a  computer  without  prior  knowledge  of  computer  languages .  Every  effort 
has  been  made  to  produce  a  system  as  machine  independent  as  possible  to 
make  implementation  on  any  large  computer  configuration  relatively  easy. 
However,  there  are  a  few  modifications  which  may  have  to  be  made. 

This  Technical  Note  provides  assistance  to  the  systems  programmer, 
with  the  task  of  implementing  OMNITAB  II,  by  pointing  out  where 
difficulties  may  occur  and  how  to  cope  with  them.  It  furthermore 
outlines  a  method  for  segmenting  the  OMNITAB  II  system  which  is  very 
large.  It  is  a  partial  documentation  of  the  OMNITAB  program. 

OMNITAB  II  is  a  large  system  requiring  a  large  computer.  Overlay 
and  segmentation  are  virtually  essential.  A  method  for  segmenting 
OMNITAB  II  is  outlined.  The  method  should  be  useful  for  many  computers. 

Keywords:  ANSI  FORTRAN,  double  precision,  general-purpose  computer 
program,  implementation  of  OMNITAB  II,  labeled  common,  machine 
independent,  OMNITAB  II,  overlay,  segmentation,  system  parameters, 
transportable  computer  programs. 

OMNITAB  II  is  a  general-purpose  interpretive  program  which  permits  direct  use  of  a 
computer  without  a  prior  knowledge  of  computer  languages.  The  OMNITAB  II  system  allows 
users  to  utilize  a  computer  to  solve  their  problems  in  an  effective  and  accurate  manner. 

OMNITAB  was  first  conceived  by  Joseph  Hilsenrath  in  the  early  1960 's  to  permit 
scientists,  engineers  and  laboratory  technicians  the  use  of  a  large  high-speed  computer  in 
the  solution  of  their  problems  without  the  assistance  of  computer  programmers.  The  system 
was  developed  and  implemented  further  by  his  co-workers,  Hilsenrath  et  al.  (1966),  for  a 
specific  computer  and  therefore  contained  subprograms  written  in  symbolic  machine  language 
and  FORTRAN.  With  the  advent  of  third  generation  computers,  it  became  necessary  to  rewrite 
the  OMNITAB  program.  Walter  J.  Gilbert,  who  undertook  this  task  in  1966,  initiated  a  number 
of  major  changes  and  new  ideas.  Since  1968  the  Statistical  Engineering  Laboratory  of  the 
Applied  Mathematics  Division  has  been  responsible  for  maintaining  and  developing  OMNITAB. 
OMNITAB  II  (Hogben  et  al.  (1970))  is  the  result  of  this  effort. 

This  paper  is  one  of  four  which  constitute  the  documentation  for  OMNITAB  II.  A  complete 
program  listing  is  given  in  Peavy  et  al.  (1970).  A  user's  guide  is  given  in  Hogben  et  al. 
(1970).  Test  problems  with  results  are  given  in  Varner  et  al.  (1970).  The  material  in  this 
note,  without  this  introduction,  is  also  stored  on  the  magnetic  tape  which  contains  the 
OMNITAB  II  master  program  and  is  available  from  National  Technical  Information  Service 
(formerly  Clearinghouse),  U.S.  Department  of  Commerce,  Springfield,  Va.  22151.  No  attempt 
is  made  here  to  carefully  define  terms,  such  as  worksheet,  which  are  unique  with  OMNITAB. 
It  is  assumed  the  reader  is  familar  with  Hogben  et  al.  (1970). 

Extensive  effort  has  been  exerted  to  produce  virtually  machine  independent  subprograms 
in  order  to  make  OMNITAB  II  transportable.  The  OMNITAB  II  program  has  been  rewritten  in  the 


American  National  Standards  Institute  (ANSI)  FORTRAN  language.  We  have  tried  to  avoid  the 
use  of  any  ANSI  FORTRAN  statement  which  can  not  be  successfully  compiled  on  any  specific 
computer.  All  output  is  in  120  character  per  line  format.  The  maximum  number  of 
alphanumeric  characters  per  variable  has  been  limited  to  three.  In  some  instances,  the 
steps  taken  have  resulted  in  a  loss  of  efficiency.  For  example,  use  of  the  UNIVAC  FORTRAN  V 
function  FID  would  have  simplified  the  programming  of  the  subroutine  RFORMT.  However, 
transportability  was  considered  far  more  important  than  machine  efficiency  and  a  few 
compromises  had  to  be  made. 

Special  attention  has  been  given  to  make  it  easy  to  change  pertinent  system  parameters 
for  implementation  on  computers  with  different  memory  size,  word  length  and  logical  input- 
output  units.  All  input-output  FORTRAN  statements  use  an  integer  variable  in  referencing  an 
I/O  unit  rather  than  a  specific  logical  unit.  These  variables  are  defined  in  one  subprogram 
(SETUP)  and  only  this  subroutine  needs  modification  if  the  present  assigned  values  are  not 
compatible  with  the  user's  particular  computer  configuration.  There  are  a  few  instances 
where  alphanumeric  characters  must  be  packed  to  the  full  capacity  of  the  computer  word 
length  (six  characters  per  word  are  assumed) .  Again,  the  systems  programmer  needs  to  modify 
only  one  subprogram  (PREPAK) ,  see  XREF3,  section  2.  These  features  have  enabled  systems 
programmers  to  successfully  implement  the  OMNITAB  II  system  on  many  different  large 
computers  (e.g.  IBM  360/50  up,  GE  625,  CDC  3800  and  6600,  Burroughs  5500  and  UNIVAC  1108) 
with  a  minimum  of  effort. 

The  OMNITAB  II  system  is  large  and,  therefore,  segmentation  is  necessary.  XREF3, 
section  3,  contains  an  outline  of  the  overlay  plan  used  with  the  National  Bureau  of 
Standards  computer.  Careful  attention  was  given  to  the  segmentation  in  order  to  minimize 
the  flip-flopping  of  segments  in  and  out  of  the  computer.  Statistics  shaving  the  most 
frequently  used  commands,  which  were  obtained  over  a  three  month  period,  were  employed  in 
determining  what  subprograms  were  to  be  in  residence  at  all  times.  Furthermore,  subprograms 
which  execute  related  instructions,  such  as  those  for  matrix  and  array  operations,  were 
grouped  together.  In  some  cases,  in  order  to  solve  the  overlay  problem,  it  was  necessary  to 
duplicate  a  subroutine  in  another  segment  under  a  different  subprogram  name.  This  was 
avoided  as  much  as  possible  and  was  kept  to  a  minimum.  Overlay  and/or  segmentation 
techniques  vary  greatly  with  different  computer  software  systems.  In  some  cases,  systems 
programmers  may  have  to  expend  considerable  effort  to  segment  OMNITAB  on  a  particular 
computer,  a  problem  which  will  be  investigated  soon.  If  the  computer  has  a  large  memory 
storage  and  operation  is  in  the  batch  mode,  overlay  may  not  be  required.  The  present  size 
of  OMNITAB  II,  assuming  a  floating-point  number  occupies  one  memory  location,  is 
approximately  90K  words. 

Needless  to  say,  such  a  large  scale  general-purpose  program  necessitates  extensive 
documentation  for  the  user  as  well  as  for  the  systems  programmer  implementing  OMNITAB  II. 
The  purpose  of  the  information  presented  in  this  publication  is  to  assist  in  the 
implementation  procedure.  The  documentation  is  in  the  form  of  FORTRAN  comment  statements. 
The  last  eight  columns  of  each  card  contain  the  subprogram  name,  e.g.  XR1,  and  the  line 
number,  e.g.  210,  as  in  Peavy  et  al.  (1970).  This  permits  easy  referencing  and  updating. 
This  outline  should  help  to  solve  most  of  the  problems  which  arise  in  implementing  OMNITAB 
II  on  a  computer.  This  documentation  was  prepared  as  a  preliminary  guide  rather  than  a 
definitive  manual.  A  systems  programmer's  manual  is  being  considered  to  complement  this 
guide.  Readers  who  experience  any  difficulties  in  the  implementation  of  OMNITAB  are  asked 
to  report  their  results  to  us  so  that  improvements  may  be  made  and  passed  on  to  others. 

OMNITAB  II,  Version  5.0  was  designed  for  use  in  the  batch  processing  mode.  However,  the 
very  nature  of  OMNITAB  makes  it  adaptable  for  remote  batch  processing  or  time-sharing. 
Walter  J.  Gilbert  adapted  an  earlier  version  of  OMNITAB  for  time-sharing  which  is  being  used 
at  the  University  of  Maryland  and  the  University  of  Rome,  Italy.  One  of  the  changes  that 
has  to  be  made  in  a  time-sharing  version  is  in  the  subroutine  ERROR,  see  Peavy  et  al. 
(1970).  The  statement  on  line  ERR  210  needs  to  be  changed  so  that  NERROR=0  instead  of 
NERR0R=NERR0R+1.  If  a  terminal  having  only  72  characters  per  line  is  used,  major  changes 
would  have  to  be  made  in  the  subprograms  which  execute  the  instructions  which  have  a 
comprehensive  automatic  printing.  See  Hogben  et  al.  (1970)  for  a  listing.  Also,  a  change 


would  have  to  be  made  in  the  command  PLOT,  although  the  command  PAGE  PLOT  could  be  modified 
and  used  to  replace  PLOT.  When  the  necessary  hardware  and  software  are  added  to  the  NBS 
computer  configuration,  a  timesharing  version  of  CMNITAB  will  be  developed. 

We  thank  David  Hogben  for  his  valuable  assistance  and  suggestions. 

The  documentation  is  grouped  as  follows : 

XREF  Page  5 

This  gives  an  outline  of  the  information  contained  in  XREF1,  XREF2,  XREF3,  XREF4,  XREF5 
and  XREF6. 

XREF1  Page  6 

A  complete  list  of  all  the  labeled  common  statements  and  the  subprograms  in  which  they 
appear  is  given.  A  list  of  subprograms  and  the  referencing  subprograms  is  provided.  Also 
included  are  the  system  library  functions  (e.g.  SIN,  COS,  etc.)  used  in  the  OMNITAB  II 
system  and  the  subprograms  which  reference  them.  Enumerated  are  the  subprograms  which 
contain  DOUBLE  PRECISION  statements  along  with  the  variables  typed  as  double  precision. 
EQUIVALENCE  statements  are  used  in  a  number  of  subprograms.  The  variables  and  their 
equivalent  variables  are  listed  along  with  the  subprograms  in  which  they  appear.  Hence,  it 
should  be  easier  to  resolve  any  difficulties  caused  by  the  use  of  (or  changes  in  the  use  of) 
DOUBLE  PRECISION  or  EQUIVALENCE  statements. 

XREF2  Page  20 

Each  CMNITAB  command  is  assigned  two  unique  values  for  internal  use.  As  the  OMNITAB 
instruction  is  scanned,  the  values  for  that  particular  command  are  stored  in  the  variables 
LI  and  L2.  These  variables  are  used  as  switches  by  the  subprogram  XECUTE  and  others  to 
control  the  flow  of  the  program  in  executing  the  instructions.  XREF2  contains  a  list  of  all 
the  command  names,  appropriate  values  for  LI  and  L2,  and  the  name  of  the  subprogram  used  to 
execute  the  command. 

XREF 3  Page  24 

The  sizes  of  the  memory  spaces  allocated  in  the  OMNITAB  II  system  for  the  worksheet  and 
scratch  areas  are  12,500  and  13,500  words  respectively.  It  may  be  necessary  to  decrease  the 
size  of  these  areas  if  a  smaller  computer  is  employed  or  segmentation  is  not  used.  In  other 
instances,  due  to  large  amounts  of  input  data,  larger  worksheet  and  scratch  areas  may  be 
needed.  Because  of  the  limitations  of  ANSI  FORTRAN,  a  number  of  changes  are  necessary  if 
the  worksheet  and  scratch  areas  are  to  be  redimensioned.  A  complete  set  of  instructions  is 
included  in  this  section  which  spell  out  the  labeled  common  statements  that  must  be 
modified,  the  subprograms  which  contain  these  labeled  common  statements,  and  a  number  of 
parameters  that  must  be  redefined  in  the  subprogram  SETUP. 

The  OMNITAB  II  program  allows  the  user  to  provide  headings  for  any  fifty  columns.  Also, 
formats  may  be  specified  for  user  controlled  input-output.  If  the  number  of  characters  per 
word  is  other  than  six,  modifications  are  necessary  before  the  FORMAT  command  can  be  used. 
Instructions  for  handling  formats  and  for  changing  the  number  of  headings  permitted  are 
given. 

An  outline  of  the  segmentation  used  at  NBS  and  the  number  of  memory  locations  needed  by 
each  segment  is  described.  There  are  fourteen  segments,  the  last  two  containing 
subsegments.  The  segments  are  labeled  as  PARTI,  PART2,  ...  ,  PART14  and  contain  the  names 
of  all  the  subprograms  and  labeled  common  areas  of  each  segment.  The  first  segment  is 
resident  in  the  computer  at  all  times  while  only  one  of  the  other  segments,  PART2  through 
PART14,  is  in  memory  storage  at  a  particular  moment.  The  size  of  each  segment  is  an 
approximation  and  dependent  on  the  fact  that  a  floating-point  number  occupies  one  memory 


location.  The  outline  given  may  not  necessarily  be  the  ideal  structure  for  overlay  on  other 
computer  configurations.  The  sizes  given  are  based  on  the  use  of  a  36  binary  bit  word. 

Four  physical  input-output  units  are  used  by  CMNITAB  II.  These  units  are  pre-assigned 
logical  units  in  the  subprogram  SETUP  which  is  called  only  once  when  CMNITAB  first  enters 
the  operating  system.  The  logical  units  may  be  reassigned,  if  necessary,  for  different 
installations. 


XREF4 


Page  34 


Bounds  are  imposed  by  CMNITAB  on  the  subprograms  evaluating  trigonometric,  exponential, 
square  root,  and  logarithmic  functions.  These  bounds  are  defined  in  DATA  statements  that 
may  need  modification  due  to  a  different  word  length  of  other  computers.  Furthermore,  the 
section  indicates  variables  which  must  be  altered  if  the  number  of  significant  digits  for  a 
real  number  is  not  eight  or  if  the  largest  integer  variable  allowed  is  less  than  2  to  the 
52nd  power. 


XREF5 


Page  37 


Four  of  the  subprograms  in  the  CMNITAB  II  system  are  BLOCK  DATA  subroutines.  This 
section  lists  the  variables  and  the  values  assigned  to  the  variables  by  the  DATA  statements 
in  the  BLOCK  DATA  subprograms.  Also  included  here,  is  the  character  set  recognized  and  used 
by  the  CMNITAB  II  system. 


XREF6 


Page  39 


As  CMNITAB  instructions  are  scanned  by  the  program,  a  thorough  check  is  made  for  errors 
and  messages  are  printed  whenever  an  error  is  detected.  Messages  are  printed  for  fatal 
errors,  arithmetic  faults  and  informative  diagnostics.  A  list  of  all  the  messages  and  their 
corresponding  number  which  is  used  when  a  subprogram  calls  the  ERROR  routine  are  included. 
See  also  Hogben  et  al.  (1970)  for  further  details. 
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XR1  3740 

XR1  3750 

XR1  3760 

XR1  3770 

XR1  3780 

XR1  3790 

XR1  3800 

XR1  3810 

XR1  3820 

XR1  3830 

XR1  3840 

XR1  3850 

XR1  3860 

XR1  3870 

XR1  3880 

XR1  3890 

XR1  3900 

XR1  3910 

XR1  39  20 

XR1  3930 

XR1  3940 

XR1  3950 

XR1  3960 

XR1  39  70 

XR1  3980 

XR1  3990 

XR1  4000 

XR1  4010 

XR1  4030 

XR1  4040 

XR1  4050 

XR1  4060 

XR1  4070 
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C  XECUTE                                              XR1  4080 

C  MOP                                                    XR1  4090 

C  XECUTE                                            XR1  4100 

C  MOVE                                                   XR1  4110 

C  MISC2  XECUTE                                      XR1  4120 

C  MPROP                                                  XR1  4130 

C  XECUTE                                            XR1  4140 

C  MRAISE                                                 XR1  4150 

C  XECUTE                                            XR1  4160 

C  MSCROW                                                 XR1  4170 

C  XECUTE                                            XR1  4180 

C  MTRIAN                                                 XR1  4190 

C  XECUTE                                            XR1  4200 

C  MTXCHK                                                 XR1  4210 

C  APRINT  ARYVEC  EXPCON  INVERT  MATRIX  MDAMAD               XR1  4220 

C  ME  I  GEN  MKRON  MMULT  MOP    MPROP  MRAISE               XR1  4230 

C  MTRIAN  MXTX   TRANSF                                XR1  4240 

C  MXTX                                                   XR1  4250 

C  XECUTE                                            XR1  4260 

C  MXTXP                                                  XR1  4270 

C  MXTX   ORTHRV                                      XR1  4280 

C  NNAME                                                  XR1  4290 

C  ASTER  OMNIT  TAPOP                                 XR1  4300 

C  NONBLA                                                 XR1  4310 

C  ASTER                                            XR1  4320 

C  NOTEPR                                                 XR1  4330 

C  OMNIT  XECUTE  XOMNIT                                XR1  4340 

C  OANOVA                                                 XR1  4350 

C  ORTHO                                            XR1  4360 

C  OCOEFF                                                 XR1  4370 

C  ORTHO                                            XR1  4380 

C  OCOVAR                                                 XR1  4390 

C  ORTHO                                            XR1  4400 

C  OMCONV                                                 XR1  4410 

C  INPUT                                            XR1  4420 

C  OMNIT                                                  XR1  4430 

C  MNITAB  OMNSYM                                      XR1  4440 

C  ONEWAY                                                 XR1  4450 

C  XECUTE                                            XR1  4460 

C  OPONE                                                  XR1  4470 

C  ORTHO                                            XR1  4480 

C  ORTHO                                                  XR1  4490 

C  XECUTE                                            XR1  4500 

C  ORTHRV                                                 XR1  4510 

C  MPROP  PROCHK                                      XR1  4520 

C  ORTPLT                                                 XR1  4530 

C  ORTHO                                            XR1  4540 

C  OUTPUT                                                 XR1  4550 

C  OMNIT                                            XR1  4560 

C  PACK                                                   XR1  4570 

C  PREPAK                                              XR1  4580 

C  PAGE                                                   XR1  4590 

C  ABRIDG  APRINT  CORREL  HISTGM  MIST   MPROP                 XR1  4600 

C  NOTEPR  OMNIT  ONEWAY  ORTHO  ORTPLT  PLOT                 XR1  4610 

C  PRINTX  RPRINT  STATIS  TWOWAY  XECUTE  XSTOP                XR1  4620 

C  PDMOTE                                                 XR1  4630 

C  XECUTE                                            XR1  4640 

C  PHYCON                                                 XR1  4670 

C  ASTER  XECUTE                                      XR1  4680 

C  PLOT                                                   XR1  4690 
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c 

XECUTE 

c 

PREPAK 

c 

ABRIDG  APRINT  HEADS  ORTHO  PRINTX  PUNCH 

c 

READX  XFORMT  XHEAD  XOMNIT 

c 

PRINTX 

c 

XECUTE 

c 

PROB 

c 

CORREL  FPPT   FPROB  OANOVA  ONEWAY 

c 

STATIS  TWOWAY 

c 

PROCHK 

c 

MPROP 

c 

PROROW 

c 

XECUTE 

c 

PUNCH 

c 

TAPOP2  XECUTE 

c 

PVTRI 

c 

MPROP 

c 

RANKO 

c 

ONEWAY 

c 

RANKS 

c 

XECUTE 

c 

RANKX 

c 

CORREL  RANKS 

c 

RCSUM 

c 

MPROP 

c 

READQ 

c 

OMNIT 

c 

READX 

c 

TAPOP2  XECUTE 

c 

REPINC 

c 

XECUTE 

c 

RESET 

c 

XECUTE 

c 

RFORMT 

c 

APRINT  HISTGM  OANOVA  OCOEFF  OCOVAR  OPONE 

c 

ORTHO  RPRINT  TWOWAY 

c 

RNDOWN 

c 

AERR   ERROR 

c 

RNJBK 

c 

FNKC 

c 

RPRINT 

c 

ABRIDG  PRINTX 

c 

SELECT 

c 

XECUTE 

c 

SET 

c 

TAPOP2  XECUTE 

c 

SETQ 

c 

OMNIT 

c 

SETUP 

c 

OMNIT 

c 

SKSYMV 

c 

PROCI  IK 

c 

SORDER 

c 

XECUTE 

c 

SORTSM 

c 

ARYVEC  MMULT  MRAISE  MT'RIAN  MXTXP  TRANSF 

c 

SPACE 

c 

XECUTE 

c 

SPINV 

c 

INVQ1K 

XR1  4700 
XR1  4710 
XR1  4720 
XR1  4730 
XR1  4740 
XR1  4750 
XR1  4760 
XR1  4770 
XR1  4780 
XR1  4790 
XR1  4800 
XR1  4810 
XR1  4820 
XR1  4830 
XR1  4840 
XR1  4850 
XR1  4860 
XR1  4870 
XR1  4880 
XR1  4890 
XR1  4900 
XR1  4910 
XR1  4920 
XR1  4930 
XR1  4940 
XR1  4950 
XR1  4960 
XR1  4970 
XR1  4980 
XR1  4990 
XR1  5000 
XR1  5010 
XR1  5020 
XR1  5030 
XR1  5040 
XR1  5050 
XR1  5060 
XR1  5070 
XR1  5080 
XR1  5090 
XR1  5100 
XR1  5110 
XR1  5120 
XR1  5130 
XR1  5140 
XR1  5150 
XR1  5160 
XR1  5170 
XR1  5180 
XR1  5190 
XR1  5200 
XR1  5210 
XR1  5220 
XR1  5230 
XR1  5240 
XR1  5250 
XR1  5260 
XR1  5270 
;XR1  5280 
XR1  5290 
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c 

STATIS 

c 

XECUTE 

c 

STMT 

c 

OMNIT 

c 

STORE 

c 

OMNIT 

c 

STORMT 

c 

MXTX 

c 

STRUVE 

c 

BESSEL 

c 

SYMV 

c 

MTRIAN  PROCHK 

c 

TAPOP 

c 

OMNIT 

c 

TAPOP2 

c 

XECUTE 

c 

THERMO 

c 

XECUTE 

c 

TPCTPT 

c 

FPPT   ONEWAY 

c 

TRANSF 

c 

MXTX 

c 

TWOWAY 

c 

XECUTE 

c 

VARCON 

c 

ASTER 

c 

VECTOR 

c 

DEFINE  ERASE  EXTREM  FUNCT 

c 

XECUTE 

c 

OMNIT 

c 

XFORMT 

c 

OMNIT 

c 

XHEAD 

c 

OMNIT 

c 

XOMNIT 

c 

OMNIT 

c 

XPND 

c 

EXPAND  REPINC 

c 

XSTOP 

c 

OMNIT  XOMNIT 

c 

c***ft**********ft**CROss  REFERENCING  OF  ! 

C 

C 

FUNCTION 

REFERENCED  BY 

C 

C 

*ALOG 

C 

FLOG 

c 

*ALOG10 

c 

FLOG 10 

c 

*ATAN 

c 

FUNCT  PROB 

c 

*COS 

c 

FCOS 

c 

*DATAN2 

c 

COMPLX 

c 

*DCOS 

c 

FDCOS 

c 

*DEXP 

c 

FDEXP 

c 

*DLOG 

XR1  5300 
XR1  5310 
XR1  5320 
XR1  5330 
XR1  5340 
XR1  5350 
XR1  5360 
XR1  5370 
XR1  5380 
XR1  5390 
XR1  5400 
XR1  5410 
XR1  5420 
XR1  5430 
XR1  5440 
XR1  5450 
XR1  5460 
XR1  5470 
XR1  5480 
XR1  5490 
XR1  5500 
XR1  5510 
XRI  5520 
XR1  5530 
XRI  5540 
XRI  5550 
XRI  5560 

MISC2  MSCROW  XRI  5570 

XRI  5580 
XRI  5590 
XRI  5600 
XRI  5610 
XRI  5620 
XRI  5630 
XRI  5640 
XRI  5650 
XRI  5660 
XRI  5670 
XRI  5680 
XRI  5690 
XRI  5700 

fSTEM  FUNCTIONS****************  XRI  5705 

XRI  5710 
XRI  5715 
XRI  5720 
XRI  5730 
XRI  5740 
XRI  5750 
XRI  5760 
XRI  5770 
XRI  5780 
XRI  5790 
XRI  5800 
XRI  5810 
XRI  5820 
XRI  5830 
XRI  5840 
XRI  5850 
XRI  5860 
XRI  5870 
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c 

FDLOG 

XR1  5880 

c 

*DSIN 

XR1  5890 

c 

FDSIN 

XR1  5900 

c 

*DSQRT 

XR1  5910 

c 

FDSQRT 

XR1  5920 

c 

*EXP 

XR1  5930 

c 

FEXP 

XR1  5940 

c 

*SIN 

XR1  5950 

c 

FSIN 

XR1  5960 

c 

*SQRT 

XR1  5970 

c 

FSQRT 

XR1  5980 

c 

*TANH 

XR1  5990 

c 

FTANH 

XR1  6000 

c 

XR1  6010 

c 

2     **********£)Qljg]j: 

PRECISIC 

DEFINED  IN 

SUBPROGRAMS**  XR1  6020 

c 

XR1  6030 

c 

SUBPROGRAM 

VARIABLE 

XR1  6040 

c 

XR1  6050 

c 

ARYVEC 

XR1  6060 

c 

SUM 

X 

XR1  6070 

c 

BEJN 

XR1  6080 

c 

A 

B 

C 

D 

E 

F 

XR1  6090 

c 

G 

P 

Q 

R 

X 

Y 

XR1  6100 

c 

Z 

XR1  6105 

c 

BESSEL 

XR1  6110 

c 

AA 

B 

BINT JO 

COMELL 

DBEJ 

DSNCOS 

XR1  6120 

c 

DXEX 

DXEXP 

E 

FDCOS 

FDEXP 

P 

XR1  6130 

c 

Q 

S 

T 

W 

X 

XEX 

XR1  6140 

c 

Y 

Z 

XR1  6145 

c 

BEZERO 

XR1  6150 

c 

A 

AA 

AB 

AC 

B 

FDSQRT 

XR1  6160 

c 

X 

Y 

XR1  6165 

c 

BEZONE 

XR1  6170 

c 

A 

B 

FDSQRT 

R 

S 

T 

XR1  6180 

c 

X 

Y 

XR1  6190 

c 

BINT JO 

XR1  6200 

c 

A 

B 

BINT JO 

C 

DBEJ 

X 

XR1  6210 

c 

Z 

XR1  6220 

c 

BJORCK 

XR1  6223 

c 

c 

D 

FDSQRT 

R 

Y 

XR1  6225 

c 

CBEI 

XR1  6230 

c 

A 

AA 

B 

C 

D 

E 

XR1  6240 

c 

F 

FDCOS 

FDEXP 

FDSIN 

FDSQRT 

G 

XR1  6250 

c 

H 

P 

Q 

R 

S 

T 

XR1  6260 

c 

U 

V 

W 

X 

Y 

Z 

XR1  6270 

c 

CBEK 

XR1  6280 

c 

A 

AA 

AB 

AC 

AD 

AE 

XR1  6290 

c 

B 

C 

D 

E 

F 

FDCOS 

XR1  6300 

c 

FDEXP 

FDLOG 

FDSIN 

FDSQRT 

G 

H 

XR1  6310 

c 

P 

Q 

R 

S 

T 

U 

XR1  6320 

c 

V 

W 

X 

Y 

Z 

XR1  6330 

c 

COMELL 

XR1  6340 

c 

A 

B 

C 

COMELL 

D 

E 

XR1  6350 

c 

FDLOG 

FDSQRT 

P 

Q 

X 

Z 

XR1  6360 

c 

COMPLX 

XR1  6370 

c 

D 

FDCOS 

FDSIN 

FDSQRT 

X 

Y 

XR1  6380 

c 

DBEJ 

XR1  6390 

c 

A 

B 

C 

D 

DBEJ 

DSNCOS 

XR1  6400 

c 

DXEXP 

E 

FDCOS 

FDEXP 

FDLOG 

FDSIN 

XR1  6410 

c 

FDSQRT 

H 

S 

T 

X 

Y 

XR1  6420 
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c 

DHRND 

c 

Z 

c 

ERRINT 

c 

AN 

BN 

CI 

CONS 

DN 

ERF 

c 

ERFC 

F 

FDEXP 

FN 

FNM1 

FNM2 

c 

FOUR 

GN 

GNM1 

GNM2 

ONE 

P 

c 

PREV 

RNBC 

SCF 

SUM 

TN 

TOLER 

c 

TRRTPI 

TWO 

ULCF 

ULCPS 

WN 

X 

c 

Y 

YSQ 

c 

FDCOS 

c 

DCOS 

DSNCOS 

DXEXP 

FDCOS 

X 

c 

FDEXP 

c 

DEXP 

DSNCOS 

DXEXP 

FDEXP 

X 

c 

FDLOG 

c 

DLOG 

FDLOG 

X 

c 

FDPCON 

c 

D 

X 

c 

FDSIN 

c 

DSIN 

DSNCOS 

DXEXP 

FDSIN 

X 

c 

FDSQRT 

c 

DSQRT 

FDSQRT 

X 

c 

FOURIA 

c 

A 

AA 

AB 

AC 

AD 

BA 

c 

BB 

FDCOS 

FDSIN 

R 

Y 

c 

FNEC 

c 

Y 

Z 

c 

GQUAD 

c 

B 

BMA 

BPA 

C 

DELGQ 

STORE1 

c 

STORE 2 

c 

LBCONS 

c 

DSNCOS 

DXEXP 

TRRTPI 

c 

MMULT 

c 

SUM 

X 

c 

MRAISE 

c 

SUM 

X 

c 

MTRIAN 

c 

SUM 

X 

c 

MXTX 

c 

AP 

c 

MXTXP 

c 

SUM 

XP 

c 

OANOVA 

c 

YSUM 

c 

OCOEFF 

c 

YSUM 

c 

ORTHO 

c 

DK2 

FDSQRT 

SUM 

YSUM 

c 

ORTHRV 

c 

XP 

c 

PROB 

c 

A 

B 

C 

FDCOS 

FDEXP 

FDLOG 

c 

FDSIN 

G 

ONE 

TA 

TB 

W 

c 

X 

c 

RFORMT 

c 

Z 

c 

SORTSM 

c 

SAVE 

SUM 

X 

c 

STRUVE 

c 

A 

B 

C 

DBEJ 

P 

Q 

c 

R 

S 

X 

Z 

XRl  6430 
XR1  6440 
XRl  6450 
XRl  6460 
XRl  6470 
XRl  6480 
XRl  6490 
XRl  6500 
XRl  6510 
XRl  6520 
XRl  6530 
XRl  6550 
XRl  6560 
XRl  6570 
XRl  6580 
XRl  6590 
XRl  6600 
XRl  6610 
XRl  6620 
XRl  6630 
XRl  6640 
XRl  6650 
XRl  6660 
XRl  6670 
XRl  6680 
XRl  6690 
XRl  6693 
XRl  6695 
XRl  6697 
XRl  6700 
XRl  6710 
XRl  6720 
XRl  6730 
XRl  6740 
XRl  6750 
XRl  6760 
XRl  6770 
XRl  6780 
XRl  6790 
XRl  6800 
XRl  6810 
XRl  6820 
XRl  6830 
XRl  6840 
XRl  6850 
XRl  6860 
XRl  6870 
XRl  6873 
XRl  6875 
XRl  6880 
XRl  6890 
XRl  6900 
XRl  6910 
XRl  6920 
XRl  6930 
XRl  6940 
XRl  6950 
XRl  6960 
XRl  6970 
XRl  6980 
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c 

THERMO 

XR1 

6990 

c 

EXDIF     EXX 

FDEXP     FDLOG 

G             QO 

XR1 

7000 

c 

Ql           Q2 

QQ         x 

XR1 

7010 

c 

TRANSF 

XR1 

7020 

c 

SUM         X 

XR1 

7030 

c 

TWOWAY 

XR1 

7040 

c 

DK2         FDSQRT 

SUM 

XR1 

7050 

c 

XR1 

7060 

c 

XR1 

7070 

c 

3         ** 

^VARIABLES  EQUIVALEK 

EACH  OTHER********************** 

XR1 

7080 

c 

XR1 

7090 

C  - 

SUBPROGRAM 

VARIABLE 

EQUIVALENCE 

VARIABLE 

EQUIVALENCE 

XR1 

7100 

c 

XR1 

7105 

c 

ALLSUB 

XR1 

7110 

c 

SCRA 

A 

Lll 

LL1 

XR1 

7120 

c 

L22 

LL2 

XR1 

7130 

c 

ARITH 

XR1 

7140 

c 

11 

11(1) 

12 

11(2) 

XR1 

7150 

c 

13 

11(3) 

14 

11(4) 

XR1 

7160 

c 

15 

11(5) 

XR1 

7170 

c 

ARYVEC 

XR1 

7180 

c 

X 

A 

XR1 

7190 

c 

BESSEL 

XR1 

7200 

c 

A(l) 

AA 

A(2001) 

B 

XR1 

7210 

c 

A(4001) 

W 

R 

RC 

XR1 

7220 

c 

IA 

IARGS 

KI 

KIND 

XR1 

7230 

c 

NR 

NRMAX 

XR1 

7240 

c 

CBEK 

XR1 

7250 

c 

SCRAT(1700) 

AA 

SCRAT(1800) 

AB 

XR1 

7260 

c 

COMPLX 

XR1 

7270 

c 

11 

IARGS (1) 

12 

IARGS (2) 

XR1 

7280 

c 

13 

IARGS  (3) 

14 

IARGS (4) 

XR1 

7290 

c 

15 

IARGS (5) 

16 

IARGS (6) 

XR1 

7300 

c 

CORREL 

XR1 

7310 

c 

A(13301) 

AVG(l) 

A(13401) 

SD(1) 

XR1 

7320 

c 

A(13401) 

T(D 

XR1 

7330 

c 

DBEJ 

XR1 

7340 

c 

SCRAT(1201) 

S 

SCRAT(1451) 

T 

XR1 

7350 

c 

FLIP 

XR1 

7360 

c 

I 

IARGS (100) 

J 

IARGS (99) 

XR1 

7370 

c 

K 

IARGS (9 8) 

KK 

IARGS  (9  7) 

XR1 

7380 

c 

M 

IARGS (96) 

MM 

IARGS (9 5) 

XR1 

7390 

c 

MMM 

IARGS (9 4) 

N 

IARGS (9 3) 

XR1 

7400 

c 

NN 

IARGS (9 2) 

A 

ARGS(l) 

XR1 

7410 

c 

FUNCT 

XR1 

7420 

c 

11 

11(1) 

12 

11(2) 

XR1 

7430 

c 

13 

H(3) 

14 

11(4) 

XR1 

7440 

c 

IFS 

XR1 

7450 

c 

11 

11(1) 

12 

11(2) 

XR1 

7460 

c 

13 

H(3) 

XR1 

7470 

c 

MMULT 

XR1 

7480 

c 

X 

A 

XR1 

7490 

c 

MRAISE 

XR1 

7500 

c 

X 

A 

XR1 

7510 

c 

MTRIAN 

XR1 

7520 

c 

X 

A 

XR1 

7530 

c 

MXTX 

XR1 

7540 

c 

A 

AP 

XR1 

7550 

c 

ONEWAY 

XR1 

7560 

c 

BLANK 

L(45) 

SLO 

L(22) 

XRI 

7570 

18 


c 

HIGH 

L(18) 

A2C1) 

A(2701) 

XR1  7580 

c 

A3  C.I) 

A(5401) 

A4(l) 

A(8101) 

XR1  7590 

c 

A5(l) 

A(10801) 

Bid) 

ACD 

XR1  7600 

c 

B2(l) 

A(541) 

B3(l) 

A(1081 

XR1  7610 

c 

B4(l) 

A(1621) 

B5(l) 

AC2161) 

XR1  7620 

c 

B6(l) 

A(2701) 

B7fl) 

A(3241) 

XR1  7630 

c 

B8(l) 

A(3781) 

B9(l) 

A(4321) 

XR1  7640 

c 

B10(l) 

A(4861) 

XR1  7650 

c 

OPONE 

XR1  7660 

c 

IIRGS(l) 

KFMT(l) 

XR1  7670 

c 

ORTHO 

XR1  7680 

c 

IIRGS(l) 

KFMT(l) 

BCD 

IB 

XR1  7690 

c 

ORTPLT 

XR1  7700 

c 

IU 

A 

XR1  7710 

c 

PLOT 

XR1  7720 

c 

TIT 

ITLE(1,6) 

TITX 

ITLE(1,5) 

XR1  7730 

c 

RC(1) 

X(l) 

PRINT 

A 

XR1  7740 

c 

XO 

XMIN 

XI 

XMAX 

XR1  7750 

c 

YO 

YMIN 

Yl 

YMAX 

XR1  7760 

c 

LHEAD 

IH 

IPR 

A(200) 

XR1  7770 

c 

PREPAK 

XR1  7780 

c 

A 

IAA 

XR1  7790 

c 

RFORMT 

XR1  7800 

c 

C(D 

LCD 

BLANK 

L(45) 

XR1  7810 

c 

PERIOD 

L(38) 

CPLUS 

L(40) 

XR1  7820 

c 

CMINUS 

L(39) 

CASTER 

L(41) 

XR1  7830 

c 

RPRINT 

XR1  7840 

c 

NWIDTH(l) 

A(1001) 

NDECS(l) 

A(1101) 

XR1  7850 

c 

NBLANK(l) 

AC1201) 

IRGS(l) 

A(1301) 

XR1  7860 

c 

NCOUNT(l) 

A(1401) 

XR1  7870 

c 

NWMAX(l) 

A(1601) 

NSIGDS(l) 

A(1701) 

XR1  7880 

c 

AL(1) 

LCD 

NF(1) 

A(1801) 

XR1  7890 

c 

NWM(l) 

AC1901) 

XR1  7900 

c 

SORTSM 

XR1  7910 

c 

X 

A 

XR1  7920 

c 

STATIS 

XR1  7930 

c 

A(101) 

ISA 

A(3226) 

SA 

XR1  7940 
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1830 

c 

TWOWAY 

24 

6 

TWOWAY 

UCHEBY 

19 

4 

ALLSUB 

XR2 

1840 

c 

UNIFOR 

RANDOM 

24 

15 

FNKC 

WRITE  TAPE 

47 

1 

TAPOP2 

XR2 

1850 

c 

MITE  TAPE  A 

47 

2 

TAPOP2 

WRITE  TAPE  B 

47 

3 

TAPOP2 

XR2 

1860 

c 

WRITE  TAPE  C 

47 

4 

TAPOP2 

WRITE  TAPE  D 

47 

5 

TAPOP2 

XR2 

1870 

c 

WRITE  TAPE  E 

47 

6 

TAPOP2 

WRITE  TAPE  F 

47 

7 

TAPOP2 

XR2 

1880 

c 

ZEROS 

BJONE 

30 

34 

BESSEL 

ZEROS  BJZERO 

30 

33 

BESSEL 

XR2 

1890 

c 

END 

XR2 
XR2 

1900 
1910 

23 


Listing  Of  XREF3 


C     XREF3 

C     VERSION  5.00     XREF3     5/15/70 

C 

C     THE  COMMENTS  CONTAINED  IN  THIS  SECTION  ARE  FOR  SYSTEM 

C     IMPLEMENTATION. 

C    l.**CHANGES  NEEDED  IF  WORK  SHEET  AND  SCRATCH  AREAS  ARE  TO  BE 

C       MODIFIED******** ************************************* *********  XR3 

C 

C     OMNITAB  HAS  A  WORK  SHEET  OF  12500  MACHINE  WORDS.   IF  THIS  IS  TO 

C     BE  INCREASED  OR  DECREASED,  THE  FOLLOWING  CORRECTIONS  WILL  BE 

C     NECESSARY. 

C 

C  LET  NSIZRC  BE  THE  SIZE  OF  WORK  SHEET  REQUIRED.  (NSIZRC 

C         MUST  BE  A  CONSTANT  AND  NOT  A  VARIABLE . ) 

C 

C  THEN  LET 

C 

C         NSIZR2  =  INTEGRAL  PART  OF  (NSIZRC  +  1000) /2 

C 

C         NSIZR4  =  INTEGRAL  PART  OF  NSIZRC/4 

C 

C         NSIZR5  =  INTEGRAL  PART  OF  (NSIZRC  +  1000) /5 

C 

C  NSIZR6  =  INTEGRAL  PART  OF  NSIZRC/5 

C 

C 

C   I  THE  STATEMENTS 

C 

C     COMMON/BLOCRC/NRC,RC(12600) 

C     EQUI VALENCE (ARGS(l) ,RC( 12 501)) 

C 

C     MUST  BE  CHANGED  TO 

C 

C     CQMMON/BLOCRC/NRC,RC(NSIZRC+100) 

C     EQUI VALENCE (ARGS(l) ,RC(NSIZRC+1)) 

C 

C     IN  THE  FOLLOWING  SUBPROGRAMS: 

C 

C     ABRIDG  ADRESS  ALLSUB  APRINT  ARITH  ARYVEC  BEGIN  BESSEL  CHANGE 

C     QIKCOL  CKIND  CMSEPA  COALES  COMPLX  CORREL  DEFINE  DIMENS  ERASE 

C     ERROR  EXCHNG  EXPAND  EXPCON  EXTREM  FIXFLO  FLIP   FPROB  FNEC 

C     FNEIC  FNKC   FRDIST  FUNCT  GENER  GQUAD  HISTGM  IFS    INTERP 

C     INVERT  ITERAT  LIST   MATRIX  MDAMAD  MEIGEN  MISC2  MKRON  MMULT 

C     MOP    MOVE   MPROP  MRAISE  MSCROW  MTRIAN  MXTX   OMNIT  ONEWAY 

C     OPONE  ORTHO  ORTPLT  OUTPUT  PDMOTE  PLOT   PRINTX  PROROW  PUNCH 

C     RANKS  READQ  READX  REPINC  RESET  RPRINT  SELECT  SET    SETQ 

C     SETUP  SORDER  SPACE  STATIS  STORE  THERMO  TRANSF  TWOWAY  VECTOR 

C     XCMNIT  XPND   XSTOP 

C 

C  II  THE  STATEMENT 

C 

C     COMMON/SCRAT/NS,NS2,A(13500) 
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XR3 

10 

XR3 

20 

XR3 

30 

XR3 

40 

XR3 

50 

XR3 

60 

XR3 

70 

XR3 

80 

XR3 

90 

XR3 

100 

XR3 

110 

XR3 

120 

XR3 

130 

XR3 

135 

XR3 

140 

XR3 

141 

XR3 

142 

XR3 

143 

XR3 

144 

XR3 

145 

XR3 

146 

XR3 

147 

XR3 

148 

XR3 

149 

XR3 

150 

XR3 

160 

XR3 

170 

XR3 

180 

XR3 

190 

XR3 

200 

XR3 

210 

XR3 

220 

XR3 

230 

XR3 

240 

XR3 

250 

XR3 

260 

XR3 

270 

XR3 

280 

XR3 

290 

XR3 

300 

XR3 

310 

XR3 

320 

XR3 

330 

XR3 

340 

XR3 

350 

XR3 

360 

XR3 

370 

XR3 

380 

XR3 

390 

XR3 

400 

XR3 

410 

XR3 

420 

C  XR3  430 

C  MUST  BE  CHANGED  TO                                       XR3  440 

C  XR3  450 

C  COMMON/SCRAT/NS,NS2,A(NSIZRC+1000)                           XR3  460 

C  XR3  470 

C  IN  THE  FOLLOWING  SUBPROGRAMS:                               XR3  480 

C  XR3  490 

C  ALLSUB  APRINT  ARYVEC  BESSEL  CBEK   CMSEPA  COALES  CCMPLX  CORREL    XR3  500 

C  DBEJ   EXPCON  FNEC   FNEIC  FNKC   FUNCT  HISTGM  INTERP  INVERT    XR3  510 

C  ITERAT  MATRIX  MDAMAD  ME  I  GEN  MISC2  MKRON  MMULT  MOP    MPROP     XR3  520 

C  MRAISE  MTRIAN  MXTX   OANOVA  OCOEFF  OCOVAR  ONEWAY  OPONE  ORTHO     XR3  530 

C  ORTPLT  PLOT   PREPAK  PROROW  RANKS  RPRINT  SELECT  SETUP  SORDER    XR3  540 

C  SORTSM  STATIS  TAPOP2  THERMO  TRANSF  TWOWAY  XSTOP                XR3  550 

C  XR3  560 

C  XR3  570 

C  THE  STATEMENT                                              XR3  580 

C  XR3  590 

C  COMMON/SCRAT/  NS,NS2,SCRAT( 13500)                            XR3  600 

C  XR3  610 

C  MUST  BE  CHANGED  TO                                        XR3  620 

C  XR3  630 

C  CCMMON/SCRAT/NS,NS2,SCRAT(NSIZRC+1000)                        XR3  640 

C  XR3  650 

C  IN  THE  FOLLOWING  SUBPROGRAMS:                               XR3  660 

C  XR3  670 

C  CBEK   DBEJ                                               XR3  680 

C  XR3  690 

C  (  THE  VARIABLE  A  WAS  CHANGED  TO  SCRAT  BECAUSE  A  WAS  USED         XR3  700 

C  IN  THE  SUBPROGRAMS)                                       XR3  710 

C  XR3  720 

C  III  THE  STATEMENTS  XR3  730 

C  XR3  735 

C  NS=13500                                                XR3  740 

C  NRC=12500                                                 XR3  750 

C  XR3  760 

C  MUST  BE  CHANGED  TO                                        XR3  770 

C  XR3  780 

C  NS=NSIZRC+1000                                           XR3  790 

C  NRC=NSIZRC                                              XR3  800 

C  XR3  810 

C  IN  THE  SUBPROGRAM  SETUP.                                  XR3  820 

C  XR3  830 

C  IV  THE  STATEMENTS  XR3  840 

C  XR3  850 

C  DIMENSION  SAC3125, 3), ISA(3125)                              XR3  860 

C  EQUIVALENCE  (A(101), ISA), (A(3226),SA),(NRMAX,NARMAX)            XR3  870 

C  XR3  880 

C  MUST  BE  CHANGED  TO                                        XR3  890 

C  XR3  900 

C  DIMENSION  SA(NSIZ4, 3), ISA(NSIZR4)                            XR3  910 

C  EQUIVALENCE  (A(101) ,ISA) , (A(NSIZR4+101) ,SA) , (NRMAX,NARMAX)       XR3  920 

C  XR3  930 

C  IN  THE  SUBPROGRAM  STATIS.                                  XR3  940 

C  XR3  950 

C   V  THE  STATEMENTS  XR3  960 

C  XR3  970 

C  NROW=201                                                XR3  980 

C  NCOL=62                                                 XR3  990 

C  XR3  1000 

C  MUST  BE  CHANGED  TO                                       XR3  1010 


25 


C  XR3  1020 

C  NROW=KROW  XR3  1030 

C  NCOL=KCOL  XR3  1040 

C  WHERE  KCOL*KROW  IS  LESS  THAN  OR  EQUAL  TO  NSIZRC  XR3  1050 

C  XR3  1060 

C  IN  THE  SUBPROGRAM  XOMNIT  XR3  1070 

C  XR3  1075 

C  VI  THE  STATEMENT  XR3  1080 

C  XR3  1090 

C  EQUIVALENCE  (A(13301) ,AVG(1)) , (A(13401) ,SD(1) ,T (1))  XR3  1100 

C  XR3  1110 

C  MUST  BE  CHANGED  TO  XR3  1120 

C  XR3  1130 

C  EQUIVALENCE  (A(NSIZRC+801) ,AVG(1)) , (A(NSIZRC+901) ,SD(1) ,T(1))  XR3  1140 

C  XR3  1150 

C  IN  THE  SUBPROGRAM  CORREL  XR3  1160 

C  XR3  1165 

CVII  THE  STATEMENT  XR3  1170 

C  XR3  1180 

C  DIMENSION  QQ(6750)  XR3  1190 

C  XR3  1200 

C  MUST  BE  CHANGED  TO  XR3  1210 

C  XR3  1220 

C  DIMENSION  QQ(NSIZR2)  XR3  1230 

C  XR3  1240 

C     IN  THE  SUBPROGRAM  THERMO  XR3  1250 

C  XR3  1260 

CVIII  THE  STATEMENTS  XR3  1270 

C  XR3  1280 

C  DIMENSION  A2 (2700), A3(2700) ,A4(2700) ,A5 (2700)  XR3  1290 

C  EQUIVALENCE  (A2(l) ,A(2701) )  XR3  1300 

C  EQUIVALENCE  (A3(l) ,A(5401)) , (A4(l) ,A(8101) ) , (A5(l) ,A(10801) )  XR3  1310 

C  DIMENSION  B1(540),B2(540),B3(540),B4(540),B5(540),B6(540),  XR3  1320 

C  1B7(540),B8(540),B9(540),B10(540)  XR3  1330 

C  EQUIVQLENCE  (B1(1),A(1)),(B2(1) ,A(541)),(B3(1),A(1081)),  XR3  1340 

C  1(B4(1),A(1621)),(B5(1),A(2161)),(B6(1),A(2701)),(B7(1),A(3241)),  XR3  1350 

C  2(B8(1),A(3781)),(B9(1),A(4321)),(B10(1),A(4861))  XR3  1360 

C  XR3  1370 

C  MUST  BE  CHANGED  TO  XR3  1380 

C  XR3  1395 

C  DIMENSION  A2(NSIZR5),A3(NSIZR5),A4(NSIZR5),A5(NSIZR5)  XR3  1400 

C  EQUIVALENCE  (A2(l) ,A(NSIZR5+1))  XR3  1410 

C     EQUIVALENCE  (A3(l) ,A(2*NSIZR5+1)) ,(A4(1) ,A(3*NSIZR5+1)) , (A5(l)  ,  XR3  1420 

C  1A(4*NSIZR5+1))  XR3  1430 
C     DIMENSION  B1(NSIZR6),B2(NSIZR6),B3(NSIZR6),B4(NSIZR6),B5(NSIZR6),  XR3  1440 

C  1B6(NSIZR6),B7(NSIZR6),B8(NSIZR6),B9(NSIZR6),B10(NSIZR6)  XR3  1450 
C     EQUIVALENCE  (Bl(l) ,A(1)  ) , (B2(l) ,A(NSIZR6+1) ) , (B3(l) ,A(2*NSIZR6+1))XR3  1460 

C  1,(B4(1),A(3*NSIZR6+1)),(B5(1),A(4*NSIZR6+1)),(B6(1),A(5*NSIZR6+1))XR3  1470 

C  2,(B7(1),A(6*NSIZR6+1)),(B8(1),A(7*NSIZR6+1)),(B9(1),A(8*NSIZR6+1))XR3  1480 

C  3,(B10(1),A(9*NSIZR6+1))  XR3  1490 

C  XR3  1500 

C     IN  THE  SUBPROGRAM  ONEWAY  XR3  1510 

C  XR3  1520 

C  2.**CHANGING  THE  NUMBER  OF  HEADINGS  PERMITTED  AND  PACKING  OF  XR3  1530 
FORMATS**** ***************************************** **********  xr3  1540 

C  XR3  1550 

2. A  INCREASING  OR  DECREASING  THE  NUMBER  OF  HEADINGS  FOR  COLUMNS.  XR3  1560 

C  XR3  1570 

2.A.1  THE  VERSION  DISTRIBUTED  ALLOWS  HEADINGS  FOR  50  COLUMNS.  ANY  XR3  1580 

50  COLUMNS  MAY  BE  HEADED.  IF  MORE  THAN  50  COLUMNS  ARE  XR3  1590 
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C  HEADED,  THE  HEADINGS  OF  THE  COLUMNS  HEADED  EARLIEST  WILL  BE  XR3  1600 

C  DELETED.  (I.E.  HEADS  ARE  STORED  IN  A  PUSH  DOWN  TABLE  WITH   XR3  1610 

C  END  OF  TABLE  PUSHED  OUT.  THIS  APPLIES  NO  MATTER  HOW  MANY   XR3  1620 

C  HEADINGS  AN  INSTALLATION  PERMITS.)                      XR3  1630 

C  XR3  1640 
C    2. A. 2  LET  NHEAD=  NUMBER  OF  COLUMN  HEADINGS  PERMITTED.  THEN  THE   XR3  1650 

C  FOLLOWING  CHANGES  HAVE  TO  BE  MADE  IN  THE  SUBPROGRAM  PREPAK.  XR3  1660 

C  XR3  1665 

C  I    DIMENSION  IFMT(12, 6), IHEAD(5, 50)                 XR3  1670 

C  XR3  1675 

C  MUST  BE  CHANGED  TO                            XR3  1680 

C  XR3  1685 

C  DIMENSION  IFMT(12, 6), IHEAD(5,NHEAD)               XR3  1690 

C  XR3  1695 

C  II   DATA  I 1/12/, LA/50/                            XR3  1700 

C  XR3  1705 

C  MUST  BE  CHANGED  TO                            XR3  1710 

C  XR3  1715 

C  DATA  II/12/,LA/NHEAD/                          XR3  1720 

C  XR3  1730 

C  (NOTE:  NHEAD  STANDS  FOR  A  INTEGER  CONSTANT)            XR3  1740 

C  XR3  1750 
C    2.B  FORMATS  MODIFICATION  FOR  COMPUTER  SYSTEMS  WHICH  DO  NOT  PACK    XR3  1760 

C       SIX  CHARACTERS  PER  MACHINE  WORD  XR3  1770 

C  XR3  1780 
C    2.B.1  FORMATS  MUST  BE  PACKED  TO  THE  MAXIMUM  NUMBER  OF  CHARACTERS  XR3  1790 

C  PER  WORD  IN  ORDER  FOR  THE  NH  CONVERSION  TO  WORK           XR3  1800 

C  XR3  1810 
C    2.B.2  LET  NCHAR=NUMBER  OF  CHARACTERS  THAT  CAN  BE  PACKED  PER  WORD  XR3  1820 

C  FOR  A  PARTICULAR  CONFIGURATION.  THEN  THE  FOLLOWING  CHANGES  XR3  1830 

C  HAVE  TO  BE  MADE  IN  THE  SUBPROGRAM  PREPAK,  WHERE  K  IS       XR3  1840 

C  DEFINED  AS.                                         XR3  1850 

C  XR3  1860 

C  K=72/NCHAR+M    M=0  IF  72/NCHAR  HAS  NO  REMAINDER  XR3  1870 

C  M=l  IF  72/NCHAR  HAS  A  REMAINDER  XR3  1880 

C  XR3  1885 

C  I    DIMENSION  IFMT(12,6),IHEAD(5,50)                 XR3  1890 

C  XR3  1895 

C  MUST  BE  CHANGED  TO                            XR3  1900 

C  XR3  1905 

C  DIMENSION  IFMT(K,  6) ,IHEAD(5,50)                 XR3  1910 

C  XR3  1915 

C  II   DATA  11/12/, LA/50/                            XR3  1920 

C  XR3  1925 

C  MUST  BE  CHANGED  TO                            XR3  1930 

C  XR3  1935 

C  DATA  II/K/,  LA/50/                            XR3  1940 

C  XR3  1945 

C  III  230  FORMAT  (12A6)                             XR3  1950 

C  XR3  1955 

C  ♦      MUST  BE  CHANGED  TO                            XR3  1960 

C  XR3  1965 

C  230  FORMAT  (KA=NCHAR)                          XR3  1970 

C  XR3  1980 

C  (NOTE:  NHEAD  AND  K  STAND  FOR  AN  INTEGER  CONSTANT)       XR3  1990 

C  XR3  2000 

C    2.B.3  BECAUSE  A  FORMAT  COMMAND  IS  TREATED  THE  WAY  ALL  OTHER    XR3  2010 

C  OMNITAB  COMMANDS  ARE  TREATED,  THE  FORMAT  COMMAND  IS      XR3  2020 

C  DUMPED  ON  THE  SCRATCH  TAPE  (OR  OTHER  I/O)  AND  REREAD     XR3  2030 

C  IN  THE  PACK  MODE  USING  THE  FORMAT  230.  THIS  IS         XR3  2040 
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C            NECESSARY  IN  ORDER  TO  MAINTAIN  MACHINE  INDEPENDENCE.  XR3  2050 

C  XR3  2060 

C  **OVERLAY  OR  SEGMENTATION***************************************  XR3  2070 

C  XR3  2080 

C    3. A    SEGMENTATION  AS  USED  ON  THE  NBS  UNIVAC  1108.  XR3  2090 

C  XR3  2100 

C          I      PART  1  IS  RESIDENT  IN  MEMORY  AT  ALL  TIMES,  XR3  2110 

C                CONTAINS  MOST  OFTEN  USED  ROUTINES,  SUBPROGRAMS  XR3  2120 

C                NEEDED  BY  MORE  THAN  ONE  OF  THE  OTHER  PARTS,  AND  XR3  2130 

C                ALL  OF  LABELED  COMMON  (EXCLUDING  LABELED  COMMON  XR3  2140 

C                FOR  BESSEL  SUBPROGRAMS.)  XR3  2150 

C  XR3  2160 

C  II     ONLY  ONE  OF  THE  OTHER  PARTS  (PART  2  THRU  PART  14)  XR3  2170 

C                IS  IN  MEMORY  AT  A  PARTICULAR  MOMENT.  XR3  2180 

C  XR3  2190 

C    3.B    OUTLINE  OF  SEGMENTATION  XR3  2200 

C  XR3  2210 

C                      PART   1  XR3  2220 

C  XR3  2230 

c  **************************************  XR3  2240 

C  *  AARGS  ACCDIG  ADRESS  AERR   ARITH  *  XR3  2250 

C  *  ASTER  BEGIN  CHKCOL  CKIND  DIMENS  *  XR3  2260 

C  *  ERROR  EXPAND  FCOS   FDCOS  FDEXP  *  XR3  2270 

C  *  FDLOG  FDPCON  FDSIN  FDSQRT  FEXP   *  XR3  2280 

C  *  FEXP2                         *  XR3  2290 

C  *  FLOG   FLOG10  FSIN   FSQRT  FTANH  **  XR3  2300 

C  *  FUNCT  GENER  HEADS  IFS    INFERR  *  XR3  2310 

C  *  INPUT  LIST   LOCATE       LOOKUP  *  XR3  2320 

C  *       MTXCHK  MXTXP  NNAME  NONBLA  *  XR3  2330 

C  *  OMCONV  OMNIT  OMNSYM  OUTPUT  PACK   *  XR3  2340 

C  *  PAGE   PHYCON  PREPAK  PROB   READQ  *  XR3  2350 

C  *  READX  REPINC  RESET  RFORMT  RNDOWN  *  XR3  2360 

C  *  RPRINT  SET    SORTSM  SPACE  SYMV   *  XR3  2370 

C  *  VARCON  VECTOR  XECUTE  XPND   XSTOP  *  XR3  2380 

C  *                              *  XR3  2390 

C  *  PLUS  LIBRARY  FUNCTIONS  (SIN, ETC)   *  XR3  2400 

C  *                              *  XR3  2410 

C  *                              *  XR3  2420 

C  *      ALL  OF  LABELED  COMMON       *  XR3  2430 

C  *                              *  XR3  2440 

C  *  ABCDEF  BLOCKA  BLOCKB  BLOCKC  BLOCKD  *  XR3  2450 

C  *  BLOCKE  BLOCKX  BLOCRC  CODE   CODETP  *  XR3  2460 

C  *  CONLB2  CONSLB  CONSTS  DCONL2  DCONLB  *  XR3  2470 

C  *  FMAT   HEADER  ICODE  ICODTP  KFMT   *  XR3  2480 

C  *  PCONST  PKSWT  QRS    SCRAT  SPRV   *  XR3  2490 

C  *  TAPE                          *  XR3  2500 

C  *                              *  XR3  2510 

C  *                              *  XR3  2520 

C  *  THERE  ARE  4  BLOCK  DATA  PROCEDURES.  *  XR3  2530 

*  THE  COMMENT  OR  2ND  CARD  HAS  THE    *  XR3  2540 

C  *  LABEL  NAME.  THE  4  PROCEDURES  ARE   *  XR3  2550 

C  *                              *  XR3  2560 

C  *  BLOCK  LBCONS  LOOKTB  PHYSIC       *  XR3  2570 

C  *                              *  XR3  2580 

Q  **************************************  y^7  2590 

C                            *  XR3  2600 

C                            *  XR3  2610 

C                            *  XR3  2620 

C                            *  XR3  2630 

C                            *  XR3  2640 
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Q  A********* A **************************************************       XR3 

C  1*1   1   1   1  1       1       *  1     *     XR3 

*  1P4RT  *  4   XR3 


CPART  *  2 
C     * 

P  ********** 

C  *  ARYVEC  * 

*  EXPCON  * 

*  MATRIX  * 

*  MDAMAD  * 
* 


*  MKRON 

*  MMULT 

*  MOP 

*  MRAISE 

*  MTRIAN 

*  MXTX 

*  STORMT 

*  TRANSF 
********** 


1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1  * 


CPART  1  5 
p  ********** 

C  *  ABRIDG  * 

*  APRINT  * 

*  FIXFLO  * 

*  NOTEPR  * 

*  PRINTX  * 

*  PUNCH  * 

*  SETQ 

*  SETUP 
STMT 
STORE 
TAPOP 
TAPOP2 

*  XFORMT 

C  *  XHEAD 

C  *  XOMNIT 
p  ********** 

C 

c 
c 
c 

CPART  *  8 

C     * 
********** 

*  COALES  * 

*  pppT   * 

*  HDIAG  * 

*  MEIGEN  * 

*  ONEWAY  * 
RANKO  * 

*  TpPTPT  * 
********** 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1*1 
* 


* 
* 
* 
* 
* 
* 
* 
* 
* 

********** 

* 


C 
C 
C 
C 
C 
C 
C 
C 

c 

C 
C 

C 


1PART  *  3 

1     * 

i  ********** 

1  *  ALLSUB  * 
1  *  CHANGE  * 
1  *  CMSEPA  *  1 


CMPARA  * 
EXCHNG  * 
FLIP   * 


PART 


1 

1  9 

1 
********** 

*  TWOWAY  * 
********** 


INTERP  *  1 
INTRP  *  1 
ITERAT  *  1 


PROROW  *  1 

SELECT  *  1 

*  SORDER  *  1 
**********  i 


PART  1  6 
********** 

*  COMPLX 

*  DEFINE 

*  ERASE 

*  EXTREM 

*  GQUAD 

*  MISC2 

*  MSCROW 

*  MOVE 

*  PDMOTE 
**********  i 

1 
1 
1 
1 
1 
1 
1 
1 


*  1 


XR3 


*  ] 
* 


*  1  **********Y]?'Z 

*  1  *  DETRNK  *XR3 

*  1  *  INVCHK  *XR3 

*  1  *  INVERT  *XR3 

*  1  *  MPROP  *XR3 

*  1  *  ORTHRV  *XR3 

*  1  *  PROCHK  *XR3 

*  1  *  PVTRI   *XR3 
*  RCSUM  *XR3 

1  *  SPINV  *XR3 

*  1  *  SKSYMV  *XR3 

*  i  **********vp2 

*  1 

*  ! 

*  1 

*  ! 

*  _. 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
********* 

* 


XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 


PART  1  7 
******** **vp3 

*  DUMMYA  *XR3 

*  DUMMYB  *XR3 

*  DUMMYC  *XR3 

*  DUMMYD  *XR3 

*  DUMMYE  *XR3 

*  DUMMYF  *XR3 
*********  *XR3 

XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
******** **XR3 

*  DHRND  *XR3 

*  ERRINT  *XR3 

*  FNEC  *XR3 

*  FNEIC  *XR3 

*  FNKC  *XR3 

*  PLOT  *XR3 

*  RNJBK  *XR3 

*  THERMO  *XR3 
**********^P3 

XR3 
XR3 


PART  *  10 
* 


2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
30  20 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 


29 


c     

C     1 

—  1 
1 

1 

1 
1 

XR3 
XR3 

3250 
3260 

CPART  1  11 

1 

PART  1  12 

1 

XR3 

3270 

C     1 

1 

1 

1 

XR3 

3280 

r    ********** 

1 

********** 

1 

XR3 

3290 

C  *  FPROB  * 

1 

*  BJORCK  * 

1 

XR3 

3300 

C  *  FRDIST  * 

1 

*  CORREL  * 

1 

XR3 

3310 

C  *  FREQCY  * 

1 

*  CSPINV  * 

1 

XR3 

3320 

C  *  HISTGM  * 

1 

*  INVCOR  * 

1 

XR3 

3330 

C  *  STATIS  * 

1 

*  MIST   * 

1 

XR3 

3340 

c   ********** 

1 

*  RANKS  * 

1 

XR3 

3350 

C 

1 

*  RANKX  * 

1 

XR3 

3360 

c 

1 

•kkkkkk-k'k'k'k 

1 

XR3 

3370 

c 

1 

1 

XR3 

3380 

c 

1 

1 

XR3 

3390 

c 

1 

1 

XR3 

3400 

c 

PART  1  13 

PART  1  14 

XR3 

3410 

c 

1 

1 

XR3 

3420 

c 

********** 

********** 

XR3 

3430 

c 

*  ABEKI 

* 

*  ORTHO  * 

XR3 

3440 

c 

*  BESSEL  * 

********** 

XR3 

3450 

c 

********** 

* 

XR3 

3460 

c 

1 

1 

XR3 

3470 

c 

1 

1 

XR3 

3480 

r             

XR3 
XR3 

3490 
3500 

C     1 

1 

l 

1 

1 

1 

1 

C     1 

1 

l 

1 

1 

1 

1 

XR3 

3510 

C     1 

1 

l 

1 

1 

1 

1 

XR3 

3520 

C     1 

1 

l 

1 

1 

1 

1 

XR3 

3530 

r   ********** 

•kkk-kk'kk'k'k'k 

********** 

1 

********** 

1 

********* 

*XR3 

3540 

C  *  BEJN   * 

*  BEZERO 

* 

*  OANOVA  * 

1 

*  OCOEFF  * 

1 

*  OCOVAR 

*XR3 

3550 

C  *  BINTJO  * 

*  BEZONE 

* 

********** 

1 

********** 

1 

********* 

*XR3 

3560 

C  *  DBEJ   * 

*  CBEI 

* 

1 

1 

XR3 

3570 

C  *  STRUVE  * 

*  CBEK 

* 

1 

1 

XR3 

3580 

c   ********** 

*  COMELL 

* 

1 

1 

XR3 

3590 

C 

*  FOURIA 

* 

1 

1 

XR3 

3600 

c 

********** 

**********      ********** 

XR3 

3610 

c 

*  OPONE  *      * 

ORTPLT  * 

XR3 

3620 

c 

**********      *^ 

:******** 

XR3 

3630 

c 

XR3 

3640 

c 

XR3 

3650 

C    3.C    APPROXIMATE  SIZE 

OF  EACH  PART  ( 

3N 

THE  NBS  COMPUTER 

XR3 

3660 

C 

THE  SIZE  AS  INDICATED  BELOW 

ASSUMES  THAT  A  MACHINE  WORD 

XR3 

3670 

C 

IS  36  BITS  LONG,  A  FLOATING 

POINT  NUMBER  USES  ONE  MEMORY 

XR3 

3680 

C 

WORD,  AND 

EACH  MACHINE  INSTRUCTION  IS  ONE  WORD 

LONG. 

XR3 

3690 

C 

XR3 

3700 

C          I 

PART 

1 

44000  WORDS 

XR3 

3710 

C 

XR3 

3720 

C 

A. 

ROUTINES      13000 

XR3 

3730 

c 

B. 

LABELED  COMMON  31000 

(THIS  INCLUDES 

A  WORK 

XR3 

3740 

c 

SHEET  OF 

1250C 

WORDS) 

XR3 

3750 

c 

C. 

NOT 

COUNTED  BUT  MUST 

BE  ADDED  ARE 

XR3 

3760 

c 

LIBRARY  FUNCTION  ] 

ROUTINES  (I.  E. 

SIN, 

COS, ETC), 

XR3 

3770 

c 

AND 

INPUT  OUTPUT 

ROUTINES.  (FOR 

THE  NBS  UNTVAC 

XR3 

3780 

c 

1108  SYSTEM  THIS 

IS 

ABOUT  4500  WORDS] 

XR3 

3790 

c 

XR3 

3800 

C         II 

PART 

2 

4100  WORDS 

XR3 

3810 

c 

XR3 

3820 

C          III     PART 

3 

4300  WORDS 

XR3 

3830 

c 

XR3 

3840 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


IV  PART  4     4000  WORDS  XR3  3850 

XR3  3860 

V  PART  5     2900  WORDS  XR3  3870 

XR3  3880 

VI  PART  6     3000  WORDS  XR3  3890 

XR3  3900 

VII  PART  7     4400  WORDS  XR3  3910 

(THIS  SEGMENT  IS  WHERE  THE  USER  MAY  XR3  3920 
ADD  HIS  SUBROUTINES)  XR3  3930 

XR3  3940 

VIII  PART  8     3700  WORDS  XR3  3950 

XR3  3960 

IX  PART  9     4000  WORDS  XR3  3970 

XR3  3980 

X  PART  10    3200  WORDS  XR3  3990 

XR3  4000 

XI  PART  11    3400  WORDS  XR3  4010 

XR3  4020 

XII  PART  12    3800  WORDS  XR3  4030 

XR3  4040 

XIII  PART  13    3300  WORDS  XR3  4050 

XR3  4060 

XIV  PART  14    4400  WORDS  XR3  4070 

XR3  4080 
4.**ASSIGNMENT  OF  INPUT  AND  OUTPUT  UNITS**************************  XR3  4090 

XR3  4100 
4. A    I/O  UNITS  FOR  OMNITAB  SYSTEM  XR3  4110 

XR3  4120 

FOUR  PHYSICAL  I/O  UNITS  ARE  USED  BY  OMNITAB.  THROUGH  OUT  ALLXR3  4130 

THE  ROUTINES  THESE  UNITS  ARE  REFERRED  BY  THE  FOLLOWING      XR3  4140 

VARIABLES  XR3  4150 

THIS  UNIT  IS  USED  FOR  READING  THE  OMNITAB  XR3  4170 

COMMANDS  AND  DATA.  (USUALLY  CARD  READER)   XR3  4180 

THE  OUTPUT  UNIT-RESULTS  AND  COMMANDS      XR3  4190 

(USUALLY  ON  LINE  PRINTER)  XR3  4200 

(62  LINES/PAGE,  120  CHAR/LINE)  XR3  4205 

3.  IPUNCH  -  THE  OUTPUT  UNIT  FOR  GENERATING  PUNCH  CARDS  XR3  4210 

4.  ISCRAT  -  A  MAGNETIC  TAPE,  DRUM,  OR  DISC,  WHICHEVER  XR3  4220 

IS  AVAILABLE  AND  FASTER.  (THIS  I/O  UNIT  IS  XR3  4230 
USED  BY  OMNITAB.)  XR3  4240 

XR3  4250 
THE  I/O  VARIABLES  ARE  ASSIGNED  LOGICAL  UNITS  IN  THE  XR3  4260 
SUBPROGRAM  CALLED  SETUP.  XR3  4270 

THEY  MAY  NEED  TO  BE  REASSIGNED  FOR  DIFFERENT  XR3  4280 

CONFIGURATIONS.  XR3  4290 

XR3  4300 
THE  PRESENT  ASSIGNMENT  IS  AS  FOLLOWS  XR3  4310 

XR3  4320 
INUNIT=5  XR3  4330 

IPRINT=6  XR3  4340 

IPUNCH=3  XR3  4350 

ISCRAT=45  XR3  4360 

XR3  4370 
4.B  PHYSICAL  TAPE  ASSIGNMENTS  FOR  OMNITAB  TAPE  COMMANDS  XR3  4380 
IN  ADDITION  USERS  MAY  REFERENCE  6  TAPE  UNITS  THROUGH  OMNITABXR3  4400 
COMMANDS  (E.G.  AS  READ  TAPE  'L' ,  WRITE  TAPE  'L' ,ETC,  WHERE  XR3  4410 
'L'  IS  EITHER  A,B,C,D,E  OR  F) .  THE  LOGICAL  UNITS  ASSIGNED  TOXR3  4420 
THESE  TYPES  ARE  7  THRU  12,  RESPECTIVELY.  IF  THESE  LOGICAL  XR3  4430 
UNITS  HAVE  TO  BE  CHANGED,  THE  FOLLOWING  CORRECTIONS  MUST  BE  XR3  4440 
MADE  IN  THE  SUBPROGRAM  TAPOP  XR3  4450 


INUNIT 


I PRINT 
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c 
c 
c 
c 
c 
c 


ITAPE=I+6 

MUST  BE  CHANGED  TO 

ITAPE=I+NEWTP  (WHERE  NEWTP  IS  THE  FIRST  LOGICAL  I/O  UNIT 
AVAILABLE  MINUS  1) 


XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
5.**DIRECTIONS  FOR  THE  USE  OF  DUMMY  A-F***************************  XR3 

XR3 

THE  OMNI TAB  COMMANDS  DUMMY  A  THROUGH  DUMMY  F  ENABLE  THE  USERXR3 

TO  EXECUTE  HIS  OWN  SUBPROGRAM(S)  THROUGH  OMNITAB.  THE  FIRST  XR3 

STATEMENT  OF  THE  SUBPROGRAM  MUST  BE  ONE  OF  THE  FOLLOWING    XR3 

SUBPROGRAM  DUMMYA 

SUBPROGRAM  DUMMYB 

SUBPROGRAM  DUMMYC 

SUBPROGRAM  DUMMYD 

SUBPROGRAM  DUMMYE 

SUBPROGRAM  DUMMYF 


IN  WRITING  THE  SUBPROGRAM(S)  THE  FOLLOWING  INFORMATION 
SHOULD  BE  FOUND  USEFUL. 


XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 


THE  FOLLOWING  IABELED  COMMON,  DIMENSION  AND  EQUIVALENCE 
STATEMENTS  SHOULD  BE  USED. 

*******************  XR3 

COMMON/BLOCRC/NRC,RC(NSIZRC+100j       **SEE  XREF3  PARTI**  XR3 

DIMENSION  ARGS(IOO)  **FOR  DEFINITION  **  XR3 

EQUIVALENCE  ARGS(l) ,RC(NSIZRC+1)       **OF  NSIZRC      **  XR3 

COMMON/SCRAT/NR,NS2,A(NSIZRC+1000)     *******************  XR3 
COMMON/BLOCKD/IARGS(100) ,KIND(100) ,ARGTAB(100) ,NRMAX, 
1  NROW,NCOL,NARGS,VWXYZ(8),NERROR 


**VMIABLE** 
RC 
NRC 
NROW 
NCOL 
NRMAX 
A 
NS 
NS2 
IARGS 

MGS 

NARGS 


XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
XR3 
ALL  INTEGER  ARGUMENTS  OF  THE  DUMMYA, . . . ,OR  XR3 
DUMMYF  INSTRUCTION  XR3 

XR3 
ALL  FLOATING  ARGUMENTS  OF  THE  DUMMYA, . . . ,  XR3 
OR  DUMMYF  INSTRUCTION  XR3 

XR3 
THE  NUMBER  OF  ARGUMENTS  OF  THE  DUMMYA, . . . ,  XR3 
OR  DUMMYF  INSTRUCTION  XR3 

XR3 


**MEANING  OF  VARIABLE** 
THE  NAME  OF  THE  WORKSHEET  AREA 
THE  SIZE  OF  THE  WORKSHEET  AREA 
THE  NUMBER  OF  ROWS  IN  THE  WORKSHEET 
THE  NUMBER  OF  COLUMNS  IN  THE  WORKSHEET 
THE  PRESENT  NO.  OF  ROWS  BEING  USED 
AVAILABLE  SCRATCH  AREA 
SIZE  OF  THE  SCRATCH  AREA 
ONE-HALF  THE  SIZE  OF  THE  SCRATCH  AREA 


4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
5090 
5100 
5110 
5120 
5130 
5150 
5160 
5170 
5180 
5190 
5200 
5230 
5240 
5250 
5260 
5270 
5280 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 
5610 
5620 
5630 
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C  KIND               KIND(I)=0  IF  ITU  ARGUMENT  IS  AN  INTEGER           XR3  5640 

C  KIND (I) =1   IF  ITH  ARGUMENT  IS  FLOATING               XR3  5650 

C  XR3  5660 

C  NERROR           NERROR=0   IF  THERE  HAVE  BEEN  NO  FATAL  ERRORS  XR3  5670 

C  NERROR=  NO.  OF  PREVIOUS  FATAL  ERRORS                 XR3  5680 

C  XR3  5690 

C  THE  FOLLOWING  STATEMENTS  SHOULD  BE  USED  FOR  CHECKING                 XR3  5700 

C  PURPOSES  AND  TO  COMPUTE  ADDRESSES  OF  COLUMNS  IN  THE                   XR3  5710 

C  WORKSHEET.                                                                                                        XR3  5720 

C  XR3  5  7  30 

C  CALL  CKIND(J)                                                                                                XR3  5740 

C  XR3  5750 

C  CKIND  IS  A  SUBPROGRAM  WHICH  CHECKS  THE  KIND  OF  ARGUMENTS  XR3  5760 

C  (INTEGER  OR  FLOATING) .J  VALUES  ARE  CHECKED.                              XR3  5770 

C  UPON  RETURN  FROM  THE  SUBPROGRAM  J  IS  SET  TO  THE  FOLLOWINGXR3  5780 

C  J=0  IF  ALL  ARGUMENTS  ARE  INTEGER                                              XR3  5790 

C  J=l  IF  ALL  ARGUMENTS  ARE  FLOATING                                            XR3  5800 

C  J=2  IF  SOME  ARE  INTEGER  AND  SOME  FLOATING                           XR3  5810 

C  XR3  5820 

C  CALL  ADRESS(I,J)                                                                                          XR3  5830 

C  XR3  5840 

C  ADRESS  IS  A  SUBPROGRAM  WHICH  COMPUTES  THE  STARTING               XR3  5850 

C               '  ADDRESS  IN  THE  WORKSHEET  OF  IARGS(I) .                                          XR3  5860 

C  UPON  RETURN  FROM  THE  SUBPROGRAM  J  IS  SET  TO  THE  FOLLOWINGXR3  5870 

C  J=-(I+NRQ  IF  ARGUMENT  IS  A  FLOATING  POINT  NUMBER     XR3  5880 

C  J=0  COLUMN  NUMBER  IS  ILLEGAL                                                      XR3  5890 

C  J=  ADDRESS  OF  A  LEGAL  COLUMN                                                      XR3  5900 

C  XR3  5910 

C  CALL  CHKCOL(J)                                                                                              XR3  5920 

C  XR3  5930 

C  CHKCOL  IS  A  SUBPROGRAM  WHICH  CHECKS  THAT  ALL  NARGS               XR3  5940 

C  ARGUMENTS  ARE  LEGAL  COLUMN  NUMBERS  AND  CONVERTS  THESE         XR3  5950 

C  COLUMN  NUMBERS  TO  THEIR  STARTING  ADRESSES  IN  THE  WORK-       XR3  5960 

C  SHEET  AND  STORES  THE  ADDRESSES  IN  IARGS(I) ,   1=1,... NARGS  XR3  5970 

C  UPON  RETURN  FROM  THE  SUBPROGRAM  J  IS  SET  TO  THE  FOLLOW INGXR3  5980 

C  J=l  IF  THERE  ARE  ANY  ILLEGAL  ARGUMENTS                                  XR3  5990 

C  J=0  IF  ALL  ARGUMENTS  ARE  LEGAL  COLUMN  NUMBERS                   XR3  6000 

C  XR3  6010 

C  CALL  ERROR(I)                                                                                                XR3  6020 

C  XR3  6030 

C  IF  THERE  ARE  ERRORS  IN  THE  USER'S  CHECKING  OF  ARGUMENTS,  XR3  6040 

C  ERROR  MESSAGES  ARE  PROVIDED  BY  CALLING  SUBPROGRAM  ERROR.  XR3  6050 

C  SEE  XREF6  FOR  EXISTING  ERROR  MESSAGES.                                        XR3  6060 

C  XR3  6070 

C  IF(NERROR.NE.O)   RETURN                                                                             XR3  6080 

C  XR3  6090 

C  THE  USER  SHOULD  NOT  PROCEED  WITH  HIS  SUBPROGRAM  IF  NERRORXR3  6100 

C  IS  NOT  ZERO,  BECAUSE  INFORMATION  IN  THE  WORKSHEET  MAY  NOTXR3  6110 

C  BE  CORRECT.                                        XR3  6120 

C  XR3  6130 

Q  ***********************WApMT|jr**********************ft******  XR3  6140 

C  XR3  6150 

C  THE  USE  OF  DUMMYA, . . .DUMMYF  IS  NOT  RECOMMENDED  FOR  NOVICEXR3  6160 

C  USERS.   THE  INFORMATION  GIVEN  ABOVE  IS  PRIMARILY  FOR  THE     XR3  6170 

C  CONVENIENCE  OF  THE  SYSTEM  PROGRAMMER.                                          XR3  6080 

C  XR3  6190 

q  ***********************************************************  XR3  6200 

C  XR3  6210 

END  XR3  6220 
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Listing  Of  XREF4 


C  XREF4  XR4  10 

C  VERSION  5.00     XREF4      5/15/70  XR4  20 

C  XR4  30 

C  l.**CHANGES  THAT  MAY  BE  NEEDED  BECAUSE  SYSTEM  FUNCTIONS  ARE  USED****  XR4  40 

C  XR4  50 

C  DIFFERENT  COMPUTERS  HAVE  DIFFERENT  SIZE  WORDS  WHICH  AFFECT  THE  XR4  60 

C  SIZE  OF  REAL  AND  DOUBLE  PRECISION  NUMBERS.  BECAUSE  OF  THIS,  AN  XR4  70 

C  UPPER  LIMIT  HAS  TO  BE  PLACED  FOR  COMPUTING  THE  EXPONENTIAL,  SINE,XR4  80 

C  AND  COSINE  FUNCTIONS.  THESE  LIMITS  MAY  HAVE  TO  BE  CHANGED  FOR  XR4  90 

C  SOME  COMPUTERS.  XR4  100 

C  XR4  110 

C  IN  ORDER  TO  MAKE  OMNITAB  AS  SYSTEM  FREE  AS  POSSIBLE,  NONE  OF  THE  XR4  120 

C  PROGRAMS  MAKE  DIRECT  REFERENCE  TO  THE  FOLLOWING  FUNCTIONS:  XR4  130 

C  ALOG,ALOG10,COS,DCOS,DEXP,DLOG,DSIN,DSQRT,EXP,SIN,SQRT,TANH,  XR4  140 

C  AND  RAISING  A  VARIABLE  TO  ANOTHER  VARIABLE.  XR4  150 

C  INSTEAD  THE  OMNITAB  SUBPROGRAMS  CALL  THE  FOLLOWING  FUNCTIONS:  XR4  160 

C  FLOG,FLOG10,FCOS,FDCOS,FDEXP,FDLOG,FDSIN,FDSQRT,FEXP,FEXP2,  XR4  170 

C  FSIN,FSQRT,FTANH,  XR4  180 

C  WHICH  CHECK  TO  SEE  IF  ARGUMENT  OF  THE  FUNCTION  IS  ILLEGAL  OR  OUT  XR4  190 

C  OF  BOUNDS.  IF  THE  FUNCTION  CAN  NOT  BE  EVALUATED  INFORMATIVE  XR4  200 

C  DIAGNOSTIC  IS  PRINTED  AND  VALUE  OF  ZERO  IS  RETURNED.  IF  THE  XR4  210 

C  ARGUMENT  IS  VALID  THE  SYSTEM  FUNCTIONS  ARE  USED.  XR4  220 

C  XR4  230 

C  THE  LIMITS  ARE  DEFINED  IN  BLOCK  DATA  SUBPROGRAM  WHICH  HAS  THE  XR4  240 

C  FOLLOWING  CARD;  XR4  250 

C  BLOCK  DATA  LBCONS  XR4  260 

C  THESE  VALUES  MAY  HAVE  TO  BE  REDEFINED  IN  SOME  CASES.  XR4  270 

C  XR4  280 

C  XTRIG  IS  NOW  SET  =  3.3E7  UPPER  ABSOLUTE  LIMIT  FOR  SIN,  COS  SINGLEXR4  290 

C  PRECISION  REAL  NUMBER  XR4  300 

C  XR4  310 

C  XEXP  IS  NOW  SET  =88.0  UPPER  LIMIT  FOR  EXP  SINGLE  PREC.  REAL  NO.  XR4  320 

C  XR4  330 

C  DSNCOS  IS  NOW  SET  =3.5D16  UPPER  ABSOLUTE  LIMIT  FOR  SIN,  COS  XR4  340 

C  DOUBLE  PRECISION  REAL  NO.  XR4  350 

C  XR4  360 

C  DXEXP  IS  NOW  SET  =  70 4. DO  UPPER  LIMIT  FOR  EXP  DOUBLE  PREC.  REAL  NOXR4  370 

C  XR4  380 

C  ER  IS  NOW  SET=l.E-8  COMPUTER  ZERO.  XR4  390 

C  XR4  400 

C  NBC  IS  NOW  SET=11  NUMBER  OF  BINARY  BITS  IN  CHARACTERISTIC  OF  A  XR4  410 

C  DOUBLE  PRECISION  NUMBER  XR4  420 

C  XR4  430 

C  NBM  IS  NOW  SET=60  NUMBER  OF  BINARY  BITS  IN  MANTISSA  OF  A  DOUBLE  XR4  440 

C  PRECISION  NUMBER  XR4  450 

C  XR4  460 

C  TRRTPI  IS  NOW  SET=1. 12837916 7095512574D0  THE  VALUE  OF  2.0/SQRT(PI)XR4  470 

C  TO  19  DECIMAL  PLACES  XR4  480 

C  XR4  490 

C  XR4  500 

C  XR4  510 

XR4  520 
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C  I  FLOG,FLOG10:  IF  X  IS  GREATER  THAN  ZERO,  LOGE  OR  LOG10  IS  COMPUTED  XR4  530 

C  IF  X  IS  ZERO  OR  NEGATIVE , INFORMATIVE  DIAGNOSTIC  IS   XR4  540 

C  PRINTED  AND  ZERO  VALUE  RETURNED                  XR4  550 

C  II  FDLOG:      SAME  AS  I  EXCEPT  RESULT  IS  DOUBLE  PRECISION  VALUE    XR4  560 

C  III  FCOS,FSIN:   IF  X  IN  ABSOLUTE  VALUE  IS  LESS  THAN  OR  EQUAL  TO  3.3E7XR4  570 

C  COS  OR  SIN  IS  EVALUATED.  OTHERWISE  AN  INFORMATIVE    XR4  580 

C  DIAGNOSTIC  IS  PRINTED  AND  ZERO  VALUE  RETURNED       XR4  590 

C  IV  FDCOS,FDSIN:  IF  X  IN  ABSOLUTE  VALUE  IS  LESS  THAN  OR  EQUAL  TO     XR4  600 

C  3.5D16,  COS  OR  SIN  IS  EVALUATED  IN  DOUBLE  PREC.      XR4  610 

C  OTHERWISE  AN  INFORMATIVE  DIAGNOSTIC  IS  PRINTED  AND  XR4  620 

C  ZERO  VALUE  RETURNED.                            XR4  630 

C  V  FEXP:       IF  X  IS  LESS  OR  EQUAL  TO  88.0  THE  EXP(X)  IS  COMPUTED  XR4  640 

C  IF  VALUE  IS  NOT  EVALUATED,  THEN  AN  INFORMATIVE      XR4  650 

C  DIAGNOSTIC  IS  PRINTED  AND  ZERO  VALUE  RETURNED       XR4  660 

C  VI  FDEXP:      EXP(X)  IS  COMPUTED  IN  DOUBLE  PREC,  IF  X  IS  NOT     XR4  670 

C  GREATER  THAN  704. DO.  OTHERWISE  ZERO  VALUE  IS  RETURNEDXR4  680 

C  AND  AN  INFORMATIVE  DIAGNOSTIC  IS  PRINTED           XR4  690 

C  VII  FEXP2:      EVALUATES  X**A.  IF  ABSOLUTE  VALUE  OF  A  IS          XR4  700 

C  GREATER  THAN  60  OR  NON  INTEGER,  FEXP2  CALLS  FLOG.   XR4  710 

CVIII  FSQRT,FDSQRT: COMPUTES  SQUARE  ROOT  OF  SINGLE  AND  DOUBLE  PREC.  XR4  720 

C  NO.  RESPECTIVELY  IF  X  IS  EQUAL  OR  GREATER  THAN  ZERO.XR4  730 

C  OTHERWISE  AN  INFORMATIVE  DIAGNOSTIC  IS  PRINTED  AND  XR4  740 

C  ZERO  VALUE  RETURNED                            XR4  750 

C  IX  FTANH:       IF  ABSOLUTE  VALUE  OF  2*X  IS  LESS  THAN  OR  EQUAL  TO  88  XR4  760 

C  THE  TANH  IS  COMPUTED.  IF  NOT,  INFORMATIVE  DIAGNOSTIC  XR4  770 

C  IS  PRINTED  AND  ZERO  RETURNED.                     XR4  780 

C  XR4  785 

C  2.**CHANGES  NEEDED  IN  BLOCK  DATA  SUBPROGRAM  WHICH  DEFINES  THE       XR4  790 

C  PHYSICAL  CONSTANTS**********************************************  XR4  800 

C  XR4  810 

C  DUE  TO  THE  DIFFERENCE  IN  THE  SIZE  OF  WORDS  IN  DIFFERENT  COMPUTERS  XR4  820 

C  THE  VALUES  OF  THE  PHYSICAL  CONSTANTS  MAY  NEED  TO  BE  CHANGED.      XR4  830 

C  THESE  VALUES  ARE  DEFINED  IN  THE  BLOCK  DATA  SUBPROGRAM  WHICH  HAS   XR4  840 

C  THE  FOLLOWING  CARD;                                       XR4  850 

C  BLOCK  DATA  PHYSIC                                        XR4  860 

C  XR4  870 

C  3.**NUMBER  OF  SIGNIFICANT  DIGITS  OF  FLOATING  POINT  NOT  EIGHT********  XR4  880 

C  XR4  890 

C  CHANGES  MUST  BE  MADE  IN  SUBPROGRAMS  ACCDIG  AND  LBCONS,  IF   XR4  900 

C  NUMBER  OF  SIGNIFICANT  DIGITS  OF  A  FLOATING  POINT  NUMBER  IS  XR4  910 

C  NOT  8.                                             XR4  920 

C  XR4  930 

C  LET  SIGDIG=  NUMBER  OF  SIGNIFICANT  DIGITS  FOR  FLOATING  POINTXR4  940 

C  NUMBER.                                 XR4  950 

C  XR4  955 

C  THEN  THE  FOLLOWING  CHANGE  MUST  BE  MADE  IN  ACCDIG          XR4  960 

C  XR4  970 

C  DATA  ADMAC/8.0/                                     XR4  975 

C  XR4  980 

C  MUST  BE  CHANGED  TO                                   XR4  985 

C  XR4  990 

C  DATA  ADMAX/SIGDIG/   (WHERE  SIGDIG  IS  A  REAL  CONSTANT)      XR4  1000 

C  XR4  1010 

C  LET  JSIGD=  NUMBER  OF  SIGNIFICANT  DIGITS  FOR  FLOATING  POINT  XR4  1020 

C  NUMBER                                   XR4  1030 

C  XR4  1035 

C  THEN  THE  FOLLOWING  CHANGE  MUST  BE  MADE  IN  LBCONS.          XR4  1040 

C  XR4  1050 

C  DATA  ISIGD/8/                                       XR4  1055 

C  XR4  1060 
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C  MUST  BE  CHANGED  TO                                   XR4  1065 

C  XR4  1070 

C  DATA  ISIGD/JSIGD/  (WHERE  JSIGD  IS  AN  INTEGER  CONSTANT)     XR4  1080 

C  XR4  1090 
C  4.**WARNING  FOR  COMPUTERS  WHERE  INTEGER  NUMBERS  MUST  BE  LESS  THAN    XR4  1100 

C  10  E9  OR  ?**32  *************************************************  XR4  1110 

C  XR4  1120 

C  CHANGES  MUST  BE  MADE  IN  SUBPROGRAMS  DHRND  AND  RPRINT,  IF    XR4  1130 

C  INTEGER  NUMBERS  ARE  LESS  THAN  2**32.                    XR4  1140 

C  XR4  1150 

C  CHANGES  INVOLVE  SCALING  INTEGERS.                       XR4  1160 

C  XR4  1170 

END  XR4  1180 
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Listing  Of  XREF5 


C     XREF5  XR5  10 

C     VERSION  5.00  XREF5     5/15/70  XR5  20 

C  XR5  30 

C  l.A  VARIABLES  DEFINED  IN  SUBPROGRAM  BLOCK**********************  XR5  40 

C  XR5  50 

C  L  IS  DEFINED  IN  LABELED  COMMON  ABCDEF.  XR5  60 

C  IT  CONTAINS  THE  ANSI  FORTRAN  CHARACTER  SET:  XR5  70 

C  THE  LETTERS  A-Z  XR5  80 

C  THE  NUMBERS  0-9  XR5  90 

C  THE  CHARACTERS  =+-*/(),.$  AND  BLANK  XR5  100 

C  THE  CHARACTER  '  IS  ALSO  INCLUDED  XR5  110 

C  XR5  120 

C  NOCARD  IS  DEFINED  IN  LABELED  COMMON  FMAT.  XR5  130 

C  IT  INITIALLY  CONTAINS  THE  WORD  OMNITAB  XR5  140 

C  XR5  150 

C  ITS  CONTENTS  ARE  CHANGED  WHEN  AN  OMNITAB  CARD  IS  XR5  160 

C  ENCOUNTERED  DURING  EXECUTION.  XR5  170 

C  XR5  180 

C  IFMTS  IS  DEFINED  IN  LABELED  COMMON  FMAT.  XR5  190 

C  IT  CONTAINS  THE  STANDARD  OMNITAB  OUTPUT  FORMAT.  XR5  200 

C  IFMTS(1)=1H(     IFMTS(2)=2H1P    IFMTS(3)=2H8E  XR5  210 

C  IFMTS(4)=3H15.   IFMTS(5)=1H6     IFMTS(6)=1H)  XR5  220 

C  XR5  230 

C  XR5  240 

C  l.B  *VARIABLES  DEFINED  IN  SUBPROGRAM  LBG3NS*********************  XR5  250 

C  XR5  260 

C  XTRIG  IS  DEFINED  IN  LABELED  COMMON  CONS LB.  XR5  270 

C  XTRIG=3.3E7  AND  IS  USED  BY  FSIN  AND  FCOS  SUB PROGRAMS. XR5  280 

C  SEE  XREF4  FOR  ADDITIONAL  NOTES.  XR5  290 

C  XR5  300 

C  XEXP  IS  DEFINED  IN  LABELED  COMMON  CONSLB.  XR5  310 

C  XEXP=88.0  AND  IS  USED  BY  FEXP  SUBPROGRAM.  XR5  320 

C  SEE  XREF4  FOR  ADDITIONAL  NOTES.  XR5  330 

C  XR5  340 

C  DSNCOS  IS  DEFINED  IN  LABELED  COMMON  DCONLB.  XR5  350 

C  DSNCOS=3.5D16  AND  IS  USED  BY  FDSIN  AND  FDCOS  SUBPROG  XR5  360 

C  SEE  XREF4  FOR  ADDITIONAL  NOTES  XR5  370 

C  XR5  380 

C  DXEXP  IS  DEFINED  IN  LABELED  COMMON  DCONLB.  XR5  390 

C  DXEXP=704.0D0  AND  IS  USED  BY  FDEXP.  XR5  400 

C  SEE  XREF4  FOR  ADDITIONAL  NOTES.  XR5  410 

C  XR5  420 

C  ER  IS  DEFINED  IN  LABELED  COMMON  CONLB2.  XR5  430 

C  ER=l.E-8  AND  IS  USED  BY  CSPINV.  XR5  440 

C  SEE  XREF4  FOR  ADDITIONAL  NOTES  XR5  450 

C  XR5  460 

C  NBC  IS  DEFINED  IN  LABELED  COMMON  DCONL2.  XR5  470 

C  NBOll  AND  IS  USED  BY  ERRINT.  XR5  480 

C  SEE  XREF4  FOR  ADDITIONAL  NOTES  XR5  490 

C  XR5  500 

C  NBM  IS  DEFINED  IN  LABELED  COMMON  DCONL2.  XR5  510 

C  NBM=60  AND  IS  USED  BY  ERRINT.  XR5  520 
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SEE  XREF4  FOR  ADDITIONAL  NOTES  XR5 

XR5 

TRRTPI  IS  DEFINED  IN  LABELED  COMMON  DCONL2.  XR5 

TRRTPI=1.128379167095512574D0  AND  IS  USED  BY  ERRINT  XR5 

SEE  XREF4  FOR  ADDITIONAL  NOTES  XR5 

XR5 

ISIGD  IS  DEFINED  IN  LABELED  COMMON  CONLB2.  XR5 

ISIGD=8  AND  IS  USED  BY  FIXFLO  AND  INFERR  XR5 

SEE  XREF4  FOR  ADDITIONAL  NOTES  XR5 

XR5 
THE  FOLLOWING  VARIABLES  ARE  DEFINED  IN  LABELED  COMMON  CONSTS.XR5 

XR5 


PI 

E 

HALFPI 

DEG 

RAD 

XALOG 


PI=3. 14159265  (THE  VALUE  OF  PI) 
E=2. 71821818  (BASE  OF  NATURAL  LOGS) 
HALFPI=1. 5707963  (VALUE  OF  PI/2) 
DEG=57. 2957795  (NUMBER  OF  DEGREES  IN  ONE  RADIAN) 
RAD=. 0174532925  (NUMBER  OF  RADIANS  IN  ONE  DEGREE) 
XALOG=38.  (EXPONENT  BOUND) 
***NOTE***  PI  AND  E  ARE  ALSO  DEFINED  IN  PHYSIC  SUBPROGRAM 


XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
XR5 
l.C     VARIABLES  DEFINED  IN  SUBPROGRAM  LOOKTB*********************  XR5 

XR5 
**THE  FOLLOWING  VARIABLES  ARE  DEFINED  BY  LABELED  COMMON  CODE**  XR5 

XR5 

IALPHA  XR5 

CONTAINS  THE  OMNITAB  NUMERIC  REPRESENTATION  FOR  THE  XR5 

CHARACTERS  A,B,C,D,E,  AND  F  XR5 

THESE  CHARACTERS  ARE  USED  BY  THE  INPUT  OUTPUT  COMMANDXR5 

XR5 

NALPHA  XR5 

CONTAINS  THE  OMNITAB  NUMERIC  REPRESENTATION  FOR  THE  XR5 

CHARACTERS  V,W,X,Y,  AND  Z  XR5 

THESE  CHARACTERS  ARE  USED  BY  THE  RESET  COMMAND       XR5 

XR5 

XR5 

l.D     *VARIABLES  DEFINED  IN  SUBPROGRAM  PHYSIC*********************  XR5 

XR5 

**THE  LABELED  COMMON  PCONST  DEFINES  THE  FOLLOWING  VARIBLES**   XR5 

P  (I)  XR5 

CONTAINS  THE  NUMERICAL  VALUES  FOR  THE  PHYSICAL      XR5 

CONSTANTS  IN  SI  AND  CGS  UNITS.  XR5 

N  (I)  XR5 

CONTAINS  THE  OMNITAB  NUMERIC  REPRESENTATION  OF  THE   XR5 

PHYSICAL  CONSTANT  NAMES  USED  IN  OMNITAB.  XR5 

XR5 

***NOTE***  SEE  SUBPROGRAM  PHYSIC  FOR  ADDITIONAL  DETAILS.      XR5 

END  XR5 


530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
925 
930 
935 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1025 
1030 
1040 
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Listing  Of  XREF6 


c 

XREF6 

XR6 

10 

c 

VERSION  5.00     XREF6     5/15/70 

XR6 

20 

c 

XR6 

30 

c 

THIS  SBC  Nov  n\)Vi\!\s  \]j  ommt.vk  ERROR  MESSAGE 

XR6 

40 

c 

XR6 

50 

c 

XR6 

60 

c 

1   A       *A*AAAAAAA*pA'pAT  ERROR  fvfESSAGES***************************** 

XR6 

70 

c 

XR6 

80 

c 

ERROR 

(1) 

XR6 

90 

c 

NAME  NOT  FOUND  IN  LIBRARY 

XR6 

100 

c 

ERROR 

(2) 

XR6 

110 

c 

ILLEGAL  STATEMENT  NUMBER 

XR6 

120 

c 

ERROR 

(3) 

XR6 

130 

c 

ILLEGAL  ARGUMENT  ON  CARD 

XR6 

140 

c 

ERROR 

(5) 

XR6 

150 

c 

COMMAND  NOT  ALLOWED  IN  REPEAT  MODE 

XR6 

160 

c 

ERROR 

(6) 

XR6 

170 

c 

STATEMENT  NUMBER  MAY  NOT  BEGIN  ANY  CARD  BETWEEN  BEGIN 

XR6 

180 

c 

AND  FINISH  CARDS 

XR6 

190 

c 

ERROR 

(7) 

XR6 

200 

c 

ILLEGAL  *STATEMENT* 

XR6 

210 

c 

ERROR 

(9) 

XR6 

240 

c 

NRMAX=0 

XR6 

250 

c 

ERROR 

(10) 

XR6 

260 

c 

I4,34H  IS  AN  ILLEGAL  NUMBER  OF  ARGUMENTS 

XR6 

270 

c 

ERROR 

(11) 

XR6 

280 

c 

COLUMN  NUMBER  TOO  BIG  OR  LESS  THAN  1 

XR6 

290 

c 

ERROR 

(12) 

XR6 

300 

c 

COMMAND  STORAGE  AREA  OVERFLOW 

XR6 

310 

c 

ERROR 

(13) 

XR6 

320 

c 

STATEMENT  NUMBER  NOT  FOUND 

XR6 

330 

c 

ERROR 

(15) 

XR6 

360 

c 

DIMENSIONED  AREA  EXCEEDS  LIMIT 

XR6 

370 

c 

ERROR 

(16) 

XR6 

380 

c 

ILLEGAL  SIZE  ROW  NUMBER 

XR6 

390 

c 

ERROR 

(17) 

XR6 

400 

c 

DEFINED  MATRIX  OVERFLOWS  WORKSHEET 

XR6 

410 

c 

ERROR 

(18) 

XR6 

420 

c 

INTEGER  ARGUMENT  LESS  THAN  -8191 

XR6 

430 

c 

ERROR 

(19) 

XR6 

440 

c 

STORED  PERFORM  STATEMENT  WILL  EXECUTE  ITSELF 

XR6 

450 

c 

ERROR 

(20) 

XR6 

460 

c 

IMPROPER  TYPE  OF  ARGUMENT 

XR6 

470 

c 

ERROR 

(21) 

XR6 

480 

c 

COMMAND  MUST  BE  STORED 

XR6 

490 

c 

ERROR 

(22) 

XR6 

500 

c 

MATRIX  IS  (NEARLY)  SINGULAR 

XR6 

510 

c 

ERROR 

(23) 

XR6 

520 

c 

INSUFFICIENT  SCRATCH  AREA 

XR6 

530 

c 

ERROR 

(24) 

XR6 

540 

c 

DEGREE  IS  LARGER  THAN  NO.  OF  NON-ZERO  WEIGHTS 

XR6 

550 

c 

ERROR 

(25) 

XR6 

560 
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C  NEGATIVE  WEIGHTS  MAY  NOT  BE  USED  XR6  570 

C  ERROR  (26)  XR6  580 

C  NUMBER  OF  COLUMNS  IS  GREATER  THAN  NUMBER  OF  ROWS  XR6  590 

C  ERROR  (27)  XR6  600 

C  FORMAT  NOT  FOUND  XR6  610 

C  ERROR  (28)  XR6  620 

C  INCORRECT  TAPE  UNIT.  COMMAND  IS  NOT  EXECUTED.  XR6  630 

C  ERROR  (29)  XR6  640 

C  NUMBER  OF  ARGUMENTS  SHOULD  BE, 1 2  XR6  650 

C  ERROR  (30)  XR6  660 

C  AN  INCREMENT  COMMAND  CAN  NOT  INCREMENT  ITSELF.  XR6  670 

C  XR6  680 
c  1#B     ***********aritHMETIC  ERROR  MESSAGES************************  XR6  690 

C  XR6  700 

C  ERROR  (101)  XR6  710 

C  NEGATIVE  ARGUMENT  TO  SORT, LOG  OR  RAISE  XR6  720 

C  ERROR  (102)  XR6  730 

C  EVALUATION  OF  EXPONENT  PRODUCES  OVERFLOW  XR6  740 

C  ERROR  (103)  XR6  750 

C  ARGUMENT  OUT  OF  BOUNDS  TO  INVERSE  FUNCTION  XR6  760 

C  ERROR  (104)  XR6  770 
C            ARGUMENT  TOO  LARGE  FOR  SIN  OR  COS,  ZERO.  RETURNED , I 4 , 6H  XR6  780 

C  TIMES  XR6  790 

C  ERROR  (105)  XR6  800 

C  BESSEL  ARGUMENTS  SCALED  TO  AVOID  OVER/UNDER  FLOW.  XR6  810 

C  RETURNED, I 4, 6H  TIMES  XR6  820 

C  ERROR  (106)  XR6  830 

C  DIVISION  BY  ZERO,  RESULT  SET=0,I4,6H  TIMES  XR6  840 

C  ERROR  (107)  XR6   850 

C  TRIG  FUNCTION  NOT  DEFINED  RESULTS  SET=0,I4,6H  TIMES  XR6  855 

C  ERROR  (108)  XR6  860 
C            ONE  OF  THE  VALUES  COMPARED  IS  ZERO,  ABSOLUTE  TOLERANCE   XR6  865 

C  WAS  USED, 14, 6H  TIMES  XR6  870 

C  ERROR  (109)  XR6  873 
C            X  FOR  ELLIPTICAL  INTEGRALS  IS  =  1.0  OR  GREATER.  RESULT  ISXR6  875 

C  SET  TO  0.0.,I4,6H  TIMES  XR6  877 

C  XR6  879 
C  i.c     *********** INpoRMATIVE  DIAGNOSTIC  MESSAGES******************  XR6   880 

C  XR6  890 

C  ERROR  (201)  XR6  900 

C  TOO  MUCH  DATA  IN  SET,  READ  OR  GENERATE,  SPILL  LOST  XR6  910 

C  ERROR  (202)  XR6  920 

C  COMMAND  NOT  ALLOWED  IN  REPEAT  MODE.  EXECUTED  BUT  NOT  XR6  930 

C  STORED  XR6  940 

C  ERROR  (203)  XR6  950 

C  VALUE  REQUESTED  IN  SHORTEN,  ACOALESCE  OR  AAVERAGE  NOT  XR6  960 

C  FOUND.  XR6  9  70 

C  ERROR  (204)  XR6  980 

C  BAD  HEAD.  COLUMN  GT  50  OR  NO  /  XR6  990 

C  ERROR  (205)  XR6  1000 

C  THIS  COMMAND  WAS  NOT  EXECUTED  BECAUSE  ITS  MEANING  WAS  XR6  1010 

C  QUESTIONABLE  XR6  1020 

C  ERROR  (206)  XR6  1030 

C  F  LESS  THAN  0 ,  SET  =  0  XR6  1040 

C  ERROR  (207)  XR6  1050 

C  NU1  OR  NU2  LESS  THAN  1  XR6  1060 

C  ERROR  (208)  XR6  1070 

C  NU1  OR  NU2  TRUNCATED  TO  INTEGER  -  XR6  1080 

C  ERROR  (209)  XR6  1090 

C  IMPROPER  TITLE  NUMBER,  ASSUMED  1  XR6  1100 

40 


c 

ERROR  ( 

c 

f-1 

ERROR  ( 
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c 
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c 
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c 

ERROR  ( 
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c 

ERROR  ( 
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c 

ERROR  ( 

c 

c 

ERROR  ( 

c 

c 

210)  XR6  1110 
NO  OF  ROWS  NOT  =  TO  COLS.  MATRIX  USED  LARGEST  SQUARE  XR6  1120 

211)  XR6  1130 
ASTERISK  STRING  IMPLYING  'THRU'  INCORRECT,  IGNORED  XR6  1140 

212)  XR6  1150 
UNNECESSARY  ARGUMENTS  IN  COMMAND  IGNORED  XR6  1160 

213)  XR6  1170 
PARTIAL  STORAGE  OF  MATRIX  XR6  1175 

214)  XR6  1180 
INSUFFICIENT  SCRATCH  AREA  XR6  1185 

215)  XR6  1190 
NRMAX  IS  NOT  LARGE  ENOUGH  TO  ALLOW  ITERATION  XR6  1200 

216)  XR6  1210 
1ST  COLUMN  OF  ISETUP  OR  ISOLATE  IS  NOT  MONOTONIC  OR  IS  XR6  1220 
CONSTANT.  XR6  1230 

217)  XR6  1240 
ITERATION  HAS  FOUND  NO  VALUES.  XR6  1250 

218)  XR6  1260 
WORKSHEET  IS  TOO  SHORT  TO  ACCOMODATE  ALL  THE  VALUES  XR6  1270 
GENERATED  BY  THIS  COMMAND.  XR6  1280 

219)  XR6  1290 
MAXMIN  HAS  FOUND  NO  EXTREMA.  XR6  1300 

220)  XR6  1310 
MAXMIN  HAS  FOUND  AND  IGNORED  A  TRIAD  OF  X'S  WITH  AT  LEASTXR6  1320 
TWO  IDENTICAL  VALUES.  XR6  1330 

221)  XR6  1340 
MORE  THAN  ONE  ARGUMENT  IN  COMMAND.  ONLY  FIRST  ONE  IS  USEDXR6  1350 

222)  XR6  1360 
FORMAT  NOT  FOUND.  READABLE  FORMAT  IS  USED  XR6  1370 

223)  XR6  1380 
ONE, SOME  OR  ALL  WEIGHTS  ARE  NEGATIVE  XR6  1390 

224)  XR6  1400 
ALL  WEIGHTS  ARE  ZERO.  COMMAND  IS  NOT  EXECUTED  XR6  1410 

225)  XR6  1420 
ARG  FOR  BESIN,BESJN,BESKN  GIVES  A  RESULT  TOO  LARGE/SMALL . XR6  1430 
COMMAND  NOT  EXECUTED.  XR6  1440 

226)  XR6  1450 
COLUMN  NOT  LONG  ENOUGH  TO  STORE  ALL  ELEMENTS.  ONLY  NROW  XR6  1460 
WILL  BE  STORED.  XR6  1470 

227)  XR6  1480 
NOT  ENOUGH  DATA  ON  COL  TO  RESTORE  MATRIX/ARRAY.  DATA  XR6  1490 
AVAILABLE  WILL  BE  USED.  XR6  1500 

228)  XR6  1510 
SUM  OF  SQRS  DO  NOT  ADD  UP-ABS.  VALUE  OF  (TOTAL -ROW- COL-  XR6  1520 
RES.) /TOTAL  EXCEEDS  5.E-7  )  XR6  1530 

229)  XR6  1540 
MORE  THAN  50  HEAD  COLUMN  COMMANDS  HAVE  BEEN  USED.  XR6  1550 

230)  XR6  1560 
ATTEMPT  TO  PROMOTE  FROM  BELOW  NRMAX.  FIRST  ARGUMENT  IS  XR6  1570 
RESET  TO  NRMAX.  XR6  1580 

231)  XR6  1590 
ATTEMPT  TO  DEMOTE  OFF  THE  WORKSHEET.  SPILL  IS  LOST.  XR6  1600 

233)  XR6  1640 
NEGATIVE  VALUE (S)  WERE  ENCOUNTERED  BY  PARTITION  FUNCTION. XR6  1650 
ZEROES  STORED.  XR6  1660 

234)  XR6  1670 
NEGATIVE  ABSOLUTE  TEMPERATURES  CONVERTED.  XR6  1680 

235)  XR6  1690 
CAUTION,  USE  EXPERIMENTALLY  ONLY.  NOT  OPTIMUM  IN  ORDER  XR6  1700 
TO  MAKE  IT  MACHINE  INDEPENDENT.  REFERENCES  -  J.B.KRUSKAL,XR6  1710 
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C  ACM, 12, 92.  AND  J.H.  HALTON,SIAM  REV., 12,1.  XR6  1720 

C  ERROR  (236)  XR6  1730 

C  COMMAND  IGNORED  -  S  BEFORE  COMMAND  NAME  MEANINGLESS  IF  XR6  1740 

C  NO  STORAGE  REQUESTED.  XR6  1750 

C  ERROR  (237)  XR6  1760 

C  NUMBER  OF  SIGNIFICANT  DIGITS  AFTER  DECIMAL  PT  HAS  BEEN  XR6  1770 

C  SET  TO, 13.  XR6  1780 

END  XR6  1790 
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